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Toward an Improved Professional Spirit Among Industrial Physicists 


HERE has been a feeling in some quarters 

lately that the American Physical Society 
does not give sufficient recognition to the in- 
dustrial physicist and even that the Society does 
not welcome papers concerned with industrial 
applications of physics. For example, it has 
been pointed out that unlike the American 
Chemical Society, the American Physical Society 
has not set up divisions devoted to the various 
aspects of physics of especial interest to indus- 
trialists. In making suggestions such as the 
formation of divisions it has, of course, been 
recognized that any attempt to improve the 
service to industrial physicists should be worked 
out through the existing organizational frame- 
work of the Society. 

Actually the Society has a very considerable 
interest in industrial physics and its Council as 
well as the Governing Board of the American 
Institute of Physics has given a great deal of 
thought to the best method of encouraging 
industrial physicists to present more papers at 
meetings and to take a more active interest in 
the Society. With the support of the American 
Physical Society, the Institute has organized 
conferences on Physics in Industry, Physics in 
the Automotive Industry, Temperature, and 
Applied Nuclear Physics, and has assisted many 
university groups in the organization of similar 
meetings. In 1937, the American Physical Society 
very generously turned over to the Institute its 
journal Physics to be developed and expanded 
into the Journal of Applied Physics as a part of 
the Institute’s program to promote industrial 
physics and other forms of applied physics. 


Thus, it is evident that the Council of the 
Society has shown a most friendly interest in 
industrial physics. The formation of separate 
divisions devoted to various subjects has been 
considered several times by the Council. The 
principal objection to separate divisions is that 
an industrial physicist, such as one in the rubber 
industry, is interested in molecular physics, 
mechanics, thermodynamics, optics, x-rays, and 
almost every other branch of physics that in- 
terests the academic physicist. In fact, it has 
been the industrialists themselves who have ob- 
jected most strenuously to this so-called “‘hori- 
zontal”’ division of the Society. 

At a meeting of the Council of the American 
Physical Society on January 23, the further 
recognition and encouragement of groups of 
physicists interested in industrial and other 
special application of physics was again dis- 
cussed. In an effort to assist these groups, a 
motion was passed authorizing the president, 
secretary, and treasurer “‘to appoint various 
prominent members of the Society, acting for the 
Council, to rouse interest in special conferences 
and group sessions at the meetings of the Ameri- 
can Physical Society, and to suggest ways and 
means whereby the American Physical Society 
might render a greater service through appropri- 
ate recognition of such groups.”’ 

Thus a mechanism has been set up by the 
Council. With the support of the industrial as 
well as academic physicists, it should go far in 
encouraging the industrial physicist and in arous- 
ing a better professional spirit among physicists. 
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The Interior of the Earth 


Viewed in Relation to Earthquake Causes 
I. What Geodosy Can Tell Us About Immediate and 
Ultimate Causes of Earthquake 


BY NICHOLAS H. HECK 
U. S. Coast and Geodetic Survey, Washington, D. C. 


INTRODUCTION 


EODESY deals with the figure of the earth, 
both as a whole and in superficial detail; 
gravity ; earthtides; wandering of the poles; isos- 
tasy; and related subjects. While books and 
articles on geodetic subjects deal to some extent 
with earthquakes, the treatment is usually inci- 
dental, and accordingly it seems appropriate that 
this series of articles should include a compre- 
hensive discussion of the relation of geodesy to 
earthquake cause. 

It is appropriate to start with the subject of 
geodesy since this can tell us only about the 
earth’s crust to a depth of about 60 miles. How- 
ever, this is a very important part of the earth’s 
interior and it is the seat of the greater part of the 
earthquakes. If more were known about the crust 
we would have the key to some of the problems 
of the deeper interior of the earth. This paper 
must deal to some extent with the surface since 
that is the place of observation. 

Earthquake cause may be separated into two 
not unrelated parts—immediate cause and ulti- 
mate cause. Geodesy has a direct and demon- 
strable relation to immediate earthquake cause, 
with conclusions based directly on geodetic meas- 
urements. The relation to ultimate cause is less 
easy to demonstrate, since, while the conclusions 
are based on measurement, their interpretation 
involves theories which are to some extent still 
in a controversial stage. 


IMMEDIATE EARTHQUAKE CAUSE 


If there were only one kind of earthquake, the 
kind which sometimes occurs in the eastern part 
of the United States and Canada, which is felt 
over very large areas but with only moderate 
intensity in the central area, and no surface 
slipping, we might still be theorizing about 
immediate earthquake cause. 
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Surface Slipping 


However, there have been numerous wide- 
spread examples of earthquakes accompanied by 
surface slipping (Fig. 1). Among these are the 
following: 


Maximum Distance over 
slipping which slip 
Earthquake Date in feet occurred 
Mino Owari, 
Japan Oct. 28, 1891 20 vert. 85 miles 
»Assam, India June 12,1897 33 vert. unknown, 
moderate 
Yakutat Bay, 
Alaska Sept. 3 and 
10, 1899 474 vert. moderate 
California April 18, 1906 =. 21 hor. 275 miles 
Imperial Val- 
ley, Calif. May 18, 1940 = 12 hor. 40 miles 


Many other examples might be cited from 
Japan, India, New Zealand, and United States, 
but this suffices to indicate that immediate 
earthquake cause is related to slipping along a 
fault-plane (Fig. 2). 








Fic. 1. San Andreas Fault at right angles to road which 
was originally continuous. Vicinity of Pt. Reyes. (Courtesy 
of G. K. Gilbert, U. S. Geological Survey.) 
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It is easy to picture slipping along the probably 
rough surface of a fault-plane as capable of 
setting up earthquake waves, but the exact 
mechanism is not obvious (Fig. 3). The problem 
can be solved only through geodetic measure- 
ments since the facts must be accurately known 
over an area on both sides of the visible fault- 
trace. Distortion of fence lines or lines of trees 
may give part of the picture but it is unlikely 
that they will be in the right place, and only 
geodetic measurements can give the needed facts. 
Since geodetic measurements are an important 
factor in this problem, a brief description of the 
method is presented. 


Geodetic Survey 


A geodetic survey is one which takes into 
account the shape of the earth. Positions of 
points marked by monuments are determined by 
triangulation and elevations of bench marks by 
precise leveling. The two processes are usually 
kept distinct, that is, the elevations of triangula- 
tion stations and the exact positions of bench 
marks are not necessarily determined. 

By means of geodetic triangulation the geo- 
graphical positions (latitudes and longitudes) of 
points marked by monuments, known as stations, 
as well as the distance and direction between 
adjacent points, are determined with great accu- 
racy (Fig. 4). The triangulation consists of a 
chain of triangles, preferably so arranged as to 
make a series of connected quadrilaterals. All 











Fic. 2. Displacement of Murchison—West Coast Road 
by movement on the White Creek Fault. The difference 
in height between the dislocated road-surfaces is 14 ft. 9 in. 
(camera pointed east). (Courtesy of New Zealand Journal 
of Science.) 





tion in this presentation. 





The Interior of the Earth 


In this issue a series of papers is presented in which several 
branches of physics, i.e., measurement, mechanics, thermodynamics, 
magnetism, and radioactivity, are applied in an effort to obtain a 
better understanding of the nature of the interior of the earth and 
the relation of this interior to earthquakes. These papers are taken 
from a symposium on this subject at the centenary celebration at 
Fordham University. The Editor wishes to express his appreciation 
to Rev. Joseph Lynch, S. J., of Fordham University, for his coopera- 
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Fic. 3. Vertical fault slip—Utah 1934. 


angles are carefully measured but only an occa- 
sional length or base is measured on the ground. 
At various points in the scheme astronomical 
observations are made at so-called Laplace sta- 
tions, which aid in keeping azimuths correct 
(Fig. 8). All triangulation in North America is 
connected and rigidly adjusted and referred to 
the so-called North American Datum. Without 
going into further detail, it suffices to state that 
no effort is spared to attain the highest possible 
accuracy (Figs. 5 and 6). This assures that if tri- 
angulation is repeated, differences in position in 
excess of small unavoidable errors are real. 

In many regions it is assumed that positions 
remain unchanged. However, in regions which 
have been mentioned where crustal movements 
have occurred, this assumption is not made. 
The best practice is to have stations close to- 
gether in the vicinity of faults where movements 
have occurred, and to repeat the triangulation 
at regular intervals, especially after an earth- 
quake in which surface movement has occurred 
(Figs. 7, 8). 

Since it is necessary to know that differences 
in position are real, the repeated triangulation 
must start outside of the region affected by the 
earthquake. Though this was not done in Cali- 
fornia in the case of the 1906 earthquake, con- 
nections have since been made in several regions 
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with the undisturbed area 
to the east. Even this 
precaution might be open 
to criticism since surface 
movements have occurred 
within the last 26 years, 
in connection with 
Nevada earthquakes of 
October 2, 1915 (Pleasant 
Valley); December 20, 
1932; January 30, 1934; 
and the Utah shock of 
March 12, 1934. 

Geodetic leveling is 
adjusted to connections 
both with mean sea level 
and with intersecting 
lines, so that in the case 
of repetitions, the deter- 
mined elevations are ac- 
curate and changes are 
real within the limits of possible precision. This 
is more frequently carried out than repeated 
triangulation. 

If it is desired to obtain the position of the 
bench marks, this can be done by a method 
known as precise traverse, which is connected 
with the triangulation. 

While geodetic observations, especially level- 
ing, were made in the case of the Assam earth- 
quake, the first such observations which could 
be used for the determination of immediate 
earthquake cause was the California earthquake 
of 1906. There was sufficient earlier triangulation 
to make it possible to determine the principal 
changes, the fault-system was relatively simple, 
the terrain was well adapted to geodetic measure- 
ment, and it is reasonable to assume that the 
observed slipping is comparable in direction and 
amount to that which occurred at the earthquake 
focus of unknown but probably moderate depth. 


Elastic Rebound Theory 


The triangulation following the earthquake did 
not extend outside the region affected by move- 
ment but this did not affect the determination 
of the immediate cause of the earthquake. The 
amount of slip was greatest at the fault trace 
and rapidly decreased with distance from it, the 
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curve of decrease being similar in form at differ- 
ent points along the trace. Harry Fielding Reid, 
Geologist and Seismologist of Johns Hopkins 
University, found that the movement can best 
be explained by his elastic rebound theory. The 
complete argument appears in The Report of 
the California Earthquake Commission! and can- 
not be given here. Briefly, the results show that, 
before the earthquake, forces were set up at the 
fault by relative displacements in opposite direc- 
tions up to 5 miles away on both sides of the fault. 
This produced an elastic strain and this in turn 
was responsible for the forces which caused the 
elastic rebound when rupture occurred. It was 
the elastic rebound which set up the earthquake 
waves. 

Very often the rebound is not exactly the right 
amount to restore equilibrium and there is tem- 
porary instability. This probably accounts for 
many aftershocks and for migration of epicenters, 
that is, occurrence of successive aftershocks in 
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Fic. 4. 104th meridian” triangulation from Colorado to 
Canadian boundary. 
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Fic. 5. Modern nine-inch direction theodolite, Parkhurst 
type. This instrument is used on first-order triangulation. 
With it the average error of closure of triangles can be 
kept to one second or less. 


different parts of the fault. It is also true that 
permanent displacement tends to decrease a 
little from that occurring at the time of the 
earthquake. 

It cannot be proved that the mechanism is the 
same for more complex fault-systems but the 
basic idea is probably approximated. 

It is also possible that what occurs at consider- 
able depths should be described in different 
terms. However for the present the interpretation 
of the meaning of changes found by geodetic 
methods must come from the elastic rebound or 
some closely related theory. 


Partial Earthquake Prediction 


A real proof of the validity of the elastic re- 
bound theory as the immediate cause of earth- 
quake would be the detection of the development 
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Fic. 6. Longitude observations at a triangulation station to determine deflection 
of vertical in prime vertical. The same instrument is used also for latitude ob- 
servations to determine deflection in meridian. 


of a state of elastic strain. This would, in effect, 
be a form of earthquake prediction, not in the 
sense of exact time and place, and maximum 
intensity, but in that of indicating for a given 
region that an earthquake is impending. Several 
successful correlations of crustal changes with 
earthquake occurrence might eventually result 
in a dependable geodetic method of such predic- 
tion. If Reid is correct, triangulation just before 
the 1906 earthquake would have shown move- 
ment of points in opposite directions on opposite 
sides of the fault-plane. This could then have 
been related to triangulation after the earth- 
quake and the amount of movement necessary 
to produce rupture could have been ascertained. 


Plan of Geodetic Observations 


Obviously the assembling of satisfactory data 
of this sort will require a long time but it is im- 
portant that there should be a program of con- 
tinuous geodetic observation. It is gratifying 
that the U. S. Coast and Geodetic Survey has 
had an annual appropriation for this purpose and 
that it is aided in planning the work by the 
Advisory Committee in Seismology of the Seis- 
mological Society of America, which recently 
succeeded the similar committee of the Carnegie 
Institution of Washington which performed this 
function for many years. 

Captain C. L. Garner, Chief, Division of 


108 


Geodesy, Coast and Geodetic Survey states? 
that: “In recent years both triangulation and 
leveling have been extended over a series of 
faults using monuments spaced at 100-foot in- 
tervals near the faults and gradually increased 
to intervals of 500 feet at distances of 5 miles. 
Traverse is also run between monuments which 
are less closely spaced. With such an arrangement 
it is possible to detect very small movements 
between adjacent monuments and to obtain 
corresponding accuracy between monuments at 
some distance apart.” 

In a number of cases levels have been rerun in 
the central regions of earthquakes in Nevada, 
Montana, and Utah. Only in the case of the Utah 
earthquake of 1934 were the previous level lines 
so placed as to give an adequate picture of the 
changes. 

While progress in the application of geodetic 
methods to the determination of the condition 
of the crust which may lead to earthquake is 
being made in a few regions, the method is in 
its infancy and it is too early to discuss results, 
except possibly in Japan. 


Complexities of Interpretation 


There is evidence from Japan that the rela- 
tions may be more complicated, especially as 
there may be block movement in place of simple 
slipping. Without attempting to define the type 
of slipping, the example given by Tsuboi,* in 
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the case of the Idu earthquake of November 26, 
1930, serves to illustrate the difficulties of inter- 
pretation. Levels were run over a line about 14 
miles long on the Idu Peninsula first in 1923-4, 
then three times in 1930, and again in 1932-33. 

There were rapid changes during 1930 es- 
pecially during a period of swarm earthquakes 
which reached 609 per day on November 25, 
then on November 26 the great earthquake 
occurred. At once the rate of change slowed 
down but it had by no means stopped in January, 
1933. No reference to later levels has been found, 
but it is evident that the deduction of conditions 
leading to earthquake is by no means simple. 
The advantage of studying this matter in Japan 
has been that the motions may be rapid, and 
geodetic measurements, especially levels, have 
become routine practice, and furthermore, the 
areas involved in an important earthquake may 
be comparatively small. 











Fic. 7. Bilby steel tower, 116 feet high, for triangulation 
observations. Inner tower supports theodolite and outer 
tower the observer and other personnel. 
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Fic. 8. Part of triangulation scheme showing base line 
connection. Heavy line is base line (measured length) and 
line with arrow is Laplace azimuth (an astronomic azimuth 
corrected for deflection of vertical). 


Non-Seismic Crustal Changes 


Surface changes resembling those of seismic 
origin may in some cases be due to some other 
cause. For example, there was a lowering of a 
nearly circular region surrounding San Jose, 
California, amounting to six feet in that city. 
The lowering is associated with changes in the 
ground-water conditions. Even in a region near 
a great fault, as in the vicinity of Taft, Cali- 
fornia, in the general region of the San Andreas 
fault, there has been shearing of oil wells and 
arching of oil pipes above the ground, unaccom- 
panied by earthquake. 


Résumé of Immediate Earthquake Cause 


It therefore appears that geodetic measure- 
ments have at least partially confirmed the idea 
and established the mechanism of elastic re- 
bound, in the form of slipping along a fault-plane. 
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Further earthquake prediction in time, place, and 
maximum intensity is far in the future, but not 
in the limited sense of determining in a place 
where surface slipping has occurred before, 
whether dangerous stresses are accumulating. 
Geodetic measurements of crustal movement 
seem to be the most promising line of attack. 
It must be realized that this is true only in limited 
regions, but there are places where an accumula- 
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Fic. 9. Tiltmeter with house. There is a small Pyrex 
saucer at bottom of pit, filled with oil. When there is tilt 
there is an angle between oil surface and cover of saucer, 
producing interference rings which give amount and 
direction of tilt. 
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tion of men and property make such deter- 
mination very valuable. 

So many unwarranted earthquake predictions 
have been sprung on the public, that the diff- 
culties and limitations of one of the most 
promising methods should not be overlooked. 


Tilt Measurement 


One method of determining tilt is by geodetic 
measurements. If different changes of elevation 
are found over an area, tilt is indicated and 
local tilt could be measured by means of leveling 
However, such methods do not lend themselves 
to continuous measurement and other means 
must be used. A tiltmeter is not a geodetic instru- 
ment except in results. It may, however, be a 
seismological instrument since certain types of 
seismograph—especially the Milne-Shaw with 
10-15 second period and fairly high sensitivity— 
make excellent tiltmeters. In Japan an instru- 
ment with a quartz fiber suspension is used. In 
the United States the method of interferometry 
has been used, with recording device so that at 
a given time the angular amount and direction 
of tilt may be known (Fig. 9). 

Interpretation of tilt measurement is compli- 
cated by the fact that there are many sources of 
tilt beside that of a tectonic nature. It seems to 
have been established in Japan that tilting is 
related to earthquake occurrence, but so far no 
positive evidence has been found in the United 
States, though there have been 83 years of con- 
tinuous observation at Berkeley, California. 


ULTIMATE EARTHQUAKE CAUSE 


We now come to the discussion of possible 
relation of geodesy to ultimate earthquake cause. 
This relation can exist only through the theory 
of isostasy, a theory which has been developed 
as the result of geodetic surveys and which will 
be outlined briefly below. The question arises 
whether the process of isostatic adjustment can 
set up stresses that may result in earthquakes. 
This question has been discussed by Reid,! 
Bowie, and Daly® but without an effort at 
exhaustive treatment. Fortunately, the particular 
form of isostatic theory adopted does not affect 
the question under discussion. 
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Definition of Isostasy 


Isostasy is a condition of equilibrium of the 
outer portion of the earth called the crust. 

According to isostasy, a condition of equal 
pressure due to the outer material of the earth’s 
crust, exists at an assumed surface at a more or 
less uniform depth below sea level, perhaps of 
the order of 60 miles. Above this surface the 
average densities of the earth’s material vary 
from place to place. According to the Pratt 
theory, densities are assumed to be average 
under the coastal planes, less than normal under 
the continents, and more than normal under the 
oceans. 


Deflection of the Vertical 


One of the early puzzles in geodesy was the 
fact that if the geographic position of a point 
(latitude and longitude) was obtained by two 
different methods, there was disagreement in 
most cases (Fig. 10). The methods referred to 
are direct astronomical observation and deriving 
the position by geodetic observations from a 
distant point whose astronomical position has 
been determined. The differences are due to 
relative deflections of the vertical. Generally the 
plumb line is deflected toward land masses and 
away from large bodies of water. These condi- 
tions are exceptionally pronounced in some areas, 
as for example, the island of Puerto Rico, where 
at two points—one on the north coast and one 
on the south coast—the plumb lines are drawn 
together 56 seconds, corresponding to a mile on 
the earth. It has been found, however, that when 
the known masses are taken into account with 
reasonable assumption of density, the computed 
deflections are in all cases greater than those 
actually observed. The distribution of densities 
called for by the theory of isostasy on the other 
hand results in generally small discrepancies 
between computed and observed deflections of 
the vertical. The chief interest in the matter on 
the part of geodesists lies in determining how 
nearly the adopted figure of the earth corre- 
sponds to the true figure. 


Gravity Anomalies 


Similarly at any point on the earth’s surface 
there is a theoretical value of gravity. This quite 
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generally differs from the observed value, but 
when corrections are applied which take into 
account the distribution of mass and density 
called for by the theory of isostasy, for most of 
the earth, the anomalies are small. There are of 
course large anomalies in certain places which re- 
main large after the isostatic reductions are 
applied. 
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Fic. 10. Relation of geoid to spheroid. Under an island 
the geoid will have an upward bulge with relation to the 
spheroid. The lines AA’ and CC’ are normal to the geoid 
and lines BB’ and DD’ are normal to the spheroid. The 
angle between each pair of normal is the deflection of the 
vertical at that point. 


Isostatic Adjustment as Earthquake Cause 


Now if the crust of the earth when first solidi- 
fied was not in isostatic adjustment, this condi- 
tion must have been reached in a comparatively 
short time. Isostasy is a condition of rest, and 
such a condition when once attained will not be 
instrumental in producing earthquakes. 

Actually there is a considerable force, though 
acting over a very long time, which tends to 
disturb the condition of isostasy—that is, the 
shifting of load on the crust through erosion and 
deposition. If isostasy is to continue there must 
be an actual transfer of material by horizontal 
flow at some depth in the opposite direction to 
the transfer at the surface. While this is plastic 
flow it is easily conceivable that stresses might 
be set up in the overlying rigid crust which might 
result in earthquakes. There is little doubt that 
a large part of the earthquake activity of the 
globe is related to mountain building, including 
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whatever goes on in the deep ocean troughs which 
are obviously parts of the adjacent mountain 
systems, whether continental or island. It is held 
by some geodesists that isostatic adjustment 
may be the cause of mountain building. In 
general, the idea is that loading or unloading at 
the surface results in vertical movement in the 
upper layers, placing the material in the lower 
layers under different conditions of temperature 
with resulting expansion or contraction. 

Now, if isostatic adjustment is the sole cause 
of mountain building, the process should be 
going on only in places where there is a change 
of load through erosion and deposition. Actually 
there are many places—especially in the Pacific 
Ocean, but also in the Atlantic and Indian 
Oceans—where there are deep troughs, earth- 
quakes, and large gravity anomalies, and where 
load transfer by erosion, etc., is either absent or 
negligible. 

Therefore, the generally accepted conclusion is 
supported that there are other causes of moun- 
tain building. These can be related to geodesy 
only in case, after disturbance of isostasy by 
such a cause, isostatic adjustment, in restoring 
equilibrium, aids the mountain building process. 


Views of Vening-Meinesz 


In this connection, Meinesz’ views should be 
considered. He points out that in the region of 
the Dutch East Indies there are, parallel to the 
main trend of the mountainous islands, deep 
troughs, a belt of negative gravity-anomalies, 
and an earthquake belt. Meinesz attributes the 
negative gravity-anomaly belt to a great lateral 
compression which buckles the crust downward 
and concentrates the light material to produce 
a greater thickness than normal. The tendency 
toward isostatic adjustment may be an important 
earthquake cause in this belt, since Visser has 
found that most of the important earthquakes 
from 1909 to 1926 lie near the negative anomaly 


belt. The relation between earthquakes and 


gravity anomalies in other regions is discussed 
by Gutenberg and Richter.® 

Cooperative gravity observations in the Andes 
of Peru and Colombia by the Coast and Geodetic 
Survey and the countries concerned resulted in 


finding several places where large gravity anom- 
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alies were associated with major earthquake 
occurrence. Further investigation is needed. 


Cause of Submarine Earthquakes 


The author noted in 1936’ that few submarine 
earthquakes occurred anywhere except near the 
deep troughs and on the ocean rises. This has 
recently been confirmed by Gutenberg after 
exhaustive investigation.’ It is reasonable to 
infer that these rises—such as the Mid-Atlantic 
Ridge—are submerged mountains. However, 
there is obviously no erosion or deposition, so 
this is another confirmation of the view that 
there must be some additional cause for mountain 
building beside transfer of load by erosion and 
deposition. 

A further argument for the presence of some 
cause or causes of mountain building other than 
isostatic adjustment is that if the part of iso- 
static adjustment that is related only to erosion 
and deposition were an important earthquake 
cause, there should be many more earthquakes 
in those portions of the earthquake belts where 
there is rapid change of load than where there 
is little or no change. This is definitely not 
the case. 


Isostatic Adjustment Without Earthquake 


It is not impossible that isostatic adjustment, 
being a very slow process, can go on without 
setting up excessive strains before a slow yield- 
ing takes place, thus causing no earthquakes. 
There is definite evidence of this in Scandinavia 
and Finland, conveniently referred to as Fen- 
noscandia. There is unmistakable evidence that 
Fennoscandia has been rising at a rate which in 
certain places has been as much as one meter 
per century.*® This is undoubtedly the result of 
isostatic adjustment from the loss of the ice load 
since the last glacial period. Earthquakes have 
been few and of ‘negligible importance except 
for a strong shock in 1303 and a series of swarm 
earthquakes numbering 136 in 1846. Gravity 
observation in the vicinity of the Mississippi 


‘delta indicates that, rapid as is the accumulation 


of load, the isostatic adjustment is keeping up 
with it. This region has been practically free 
from earthquake in historic times, the only one 


recorded being a moderate one on October 19, 
1930. 
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Isostatic Adjustment as the Cause of 
Special Earthquakes 


It would seem likely that a similar adjustment 
to that in Fennoscandia should be going on in 
North America in the general region of the Great 
Lakes. There is some evidence of slow uplift.’ 
However, in this case there are some earthquakes 
of a rather unusual type, considering the earth 
as a whole. The earthquake of November 1, 1935, 
centered north of Lake Huron at Timiskaming, 
on the Canadian shield. There have been other 
widespread earthquakes along the Saint Lawrence 
Valley and, though there appears to be a line of 
weakness, the earthquakes may well be associ- 
ated with isostatic adjustment. It is possible that 
the crust has exceptional strength and does not 
yield readily to long-continued strain. 


Conclusions 


The following conclusions may then be drawn 
regarding the relation of isostatic compensation 
and ultimate earthquake cause: 

1. If isostatic compensation processes follow 
transfer of load by erosion and deposition and 
are set up in no other way, they cannot be a 
principal cause of earthquakes. 

2. With other causes acting to disturb the 
condition of isostasy, isostatic adjustment, in 
restoring equilibrium, may aid the process of 
mountain building and, therefore, may be related 
to earthquake cause. | 

3. Important isostatic adjustments may take 
place practically without earthquakes, as in 
Fennoscandia. 

4. In 


structed, possibly because the crust is too rigid 


other cases the process may be ob- 
to yield until the point of rupture is reached, 
and this may be the cause of widespread earth- 
quakes with moderate epicentral energy, as in 
eastern North America. 

5. The exact relation of gravity anomalies and 
earthquakes is not clear but it is certain that 
earthquakes occur in definite relation to belts 
of negative anomalies associated with deep ocean 
troughs. 

The discussion of continental changes indicates 
that the verification of their reality depends 
largely on geodetic observations. Some conclu- 


VOLUME 14, MARCH, 1943 





sions regarding North America have been drawn 
from levels of the Great Lakes,® from tidal 
observations along the coast, and from tilt ob- 
servations by J. P. Delaney, S.J., at Canisius 
College, Buffalo, New York, but these conclu- 
sions are open to criticism unless precise geodetic 
measurements are made. To be entirely satis- 
factory, levels of continental extent should be 
completed in a comparatively short time with- 
out adjustment to former level nets and repeated 
at intervals. The cost would perhaps be consider- 
able but my purpose is to point out possibilities 
of geodetic methods rather than to outline and 
recommend a program. 

The deep-focus earthquake has greatly dis- 
turbed the isostasist who, it seems, will have to 
revise his view that pressure is uniform below 
the depth of compensation. However, in dis- 
cussing the matter, isostasists do not seem to 
me to pay sufficient attention to the fact that 
such earthquakes occur only in limited regions, 
the belts are much smaller than those of other 
earthquakes, and they appear not to occur at all 
elsewhere. We cannot be too positive because 
the accurate computations of depth have not 
been carried out long. Generalities about the 
condition of the crust as a whole are not suited 
to this case. It should also be pointed out that 
while isostatic adjustment is held to be world- 
wide, earthquakes are far from universal. 

Daly in his book, Strength and Structure of the 
Earth,> has raised other difficulties in connection 
with the details of isostatic adjustment but dis- 
cussion of these is outside the scope of this 
paper. 

It seems likely that earth-tides and polar 
shifts can be earthquake causes, if at all, only 
in the sense of trigger forces. Efforts have been 
made to find a more direct relation but without 
result. 

The place of geodesy with reference to the 
subject of this series is then to contribute infor- 
mation about the movements of the crust rather 
than to tell about its composition, though gravity 
determinations and the theory of isostasy give 
valuable indications of the distribution of densi- 
ties. It furnishes basic information which cannot 
be neglected in the eventual solution of the 
problems.’ 
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Discussion of “The Interior of the Earth Viewed in Relation to 
Earthquakes Causes. I. What Geodosy Can Tell Us About 
Immediate and Ultimate Causes of Earthquake” 


By M. KING HUBBERT 


N independent method of analysis which 

corroborates Captain Heck’s conclusion 
that isostasy cannot be a primary cause of 
mountain building is afforded by considering the 
energy transformations involved. 

By the isostatic hypothesis of mountain 
building, as an existing mountain is eroded down 
and its mass redistributed, the achievement of a 
new equilibrium is supposed to produce a new 
mountain. 

A mountain constitutes a storage battery of 
potential energy in the earth’s gravitational 
field. Erosion, on the other hand, is a dissipative 
process, and as successive elements of the mass 
of the mountain are moved downhill by running 
water and other agents of erosion, this potential 
energy is almost completely converted into heat 
at the temperature of the surroundings, and 
eventually leaves the earth as long wave-length 
radiation. This constitutes a thermodynamically 
irreversible transformation, and as the mountain 
is reduced to a peneplane its excess of potential 
energy is almost completely dissipated and is 
not avatlable as a supply of potential energy for a 
new mountain. 

It is physically impossible to redistribute the 
mass of an existing mountain by the mechanism 
of erosion in such a manner as to produce a 
new mountain. In fact to do so would be equiva- 
lent to perpetual motion. Once a mountain is 
eroded to a peneplane, its “‘battery,”’ so to speak, 
is completely discharged and if any other 
mountain is subsequently to be built it must be 
from an entirely new supply of energy. 

Another related fact is the comparative time 
scales of isostatic adjustment and the pene- 
planation of a mountain. Both processes repre- 
sent negative exponential approaches to a final 
state of the form 


y=yo exp [—A(t—to) ]. (1) 


where yo is the initial displacement at time fo of 
the process from its final state, y the displace- 
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ment at time /, and \ the exponential constant, 
differing in value for each of the two processes. 
This can be expressed more informatively if 
we utilize the half-life period T of the process, 
defined as the amount of time required for y to 
become equal to yo/2. From this definition 


yo/2 = yee? (2) 
from which 
log. 2 0.693 


\= ——, (3) 
YF? 


Substituting the value of \ into Eq. (1) then 


gives 
t—lo 
y=Yo exp | -0.693(—) | (4) 


For isostatic adjustments the value of T as 
obtained from the post-glacial uplifts is of the 
order of magnitude of 10‘ years. For the pene- 
planation of a mountain system the half-life 
period is not too well known, but is of the order 
of 10’ years. Using these figures as approxima- 
tions, we then see that the time required to 
achieve any degree of isostatic equilibrium is of 
the order of 10-* of that required for the corre- 
sponding approach to erosional completion. 

The result of this is that a mountain under- 
going erosion should be continuously in isostatic 
equilibrium throughout the process if undis- 
turbed by other non-isostatic causes. Conversely, 
any disturbance to such a system other than a 
simple exponential approach to peneplanation 
must necessarily arise from non-isostatic causes. 

The fact that such disturbances have occurred 
repeatedly during the two billion or more years 
of geological history indicates positively that 
some large supply of energy—probably deep- 
seated within the earth—is being dissipated. 
This is the major fact of terrestrial diastrophism 
compared with which isostatic adjustment is but 
a minor and relatively unimportant phenomenon. 
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The Interior of the Earth 
Viewed in Relation to Earthquake Causes 


Il. What Radieactivity Tells Us About the Interior of the Earth 
By VicTOR F. HEss 


Fordham University, New York, New York 


T is a well-established fact that all constituents 
of the earth’s crust contain minute quantities 

of uranium, radium, thorium, and the other 
radioactive elements. A quantitative knowledge 
of their distribution among the typical primary 
and secondary rocks is of great importance for 
the geologist as well as for the seismologist. 
The methods of quantitative determination of 
the radioactive matter in rocks have been much 
improved within the last decade and we are now 
in a better position to study the role of the 
radioactive elements in the thermal history of 
our planet than ever before. This is actually the 
most important problem. We know that radium, 
thorium, and all other members of the three 
families of radioactive elements produce heat 
during their decay. Most of this heat (about 
85 percent) is contributed by the a-rays, the 
rest by 6- and y-rays, and is developed in the 
matter in which these rays are absorbed. 

This constant emission of energy in all rocks 
necessitated a radical revision of the former 
theory of the cooling of our planet. Up to about 
1930 most authors were unanimous in expressing 
the view that the’earth would steadily increase 
its temperature within geological times if the 
amount of heat calculated from the amounts of 
radioactive matter found in the surface rocks 
were produced throughout the mantle of the 
earth. The core, on the other hand, considered 
as consisting mostly of iron and nickel, was 
“believed to contain much less, if any, radioactive 
matter, since iron meteorites were always found 
to be almost free from radium. 

If we assume that the earth is in a state of 
thermal equilibrium and if the older determina- 
tions of the radioactivity of rocks by R. J. 
Strutt, J. Joly, and others, are considered as 
correct, one has to assume that the radioactive 
substances in the rocks are confined to a super- 
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ficial layer of about ten miles in thickness. 
Arthur Holmes! introduced the hypothesis that 
the concentration of radioactive matter in rocks 
diminishes with depth according to an expo- 
nential law. 

There is no direct experimental evidence 
supporting this view. The deepest bore holes, 
reaching to a depth of approximately four 
kilometers, are naturally insufficient to give 
unambiguous evidence. G. Trovato? reported 
that in the lava of Etna, a volcano in Sicily, 
the radioactivity of fresh lava deposits was 
smaller* than that of old lava. Samples taken 
from the deepest bore holes do not support the 
view that radioactivity diminishes with depth. 
On the other hand, it must be kept in mind 
that even bore holes of four and more kilo- 
meters depth are mere pin pricks considering 
the dimensions of the crust and mantle of our 
globe. 

Another possibility was to assume that when 
radioactive elements are deeply buried in the 
earth, their decay and their heat production are 
slowed up or completely stopped at sufficiently 
high pressures and temperatures. 

Here again we have no experimental evidence 
supporting this view. Temperature changes 
between about — 200°C and +2500°C, pressures 
of enormous magnitudes, and the most powerful 
magnetic fields do not, as far as we know, 
influence the rate of disintegration and the heat 
production of the known radioactive elements. 
It is also well known that the rate of radioactive 
decay of any given radioactive element does not 
change with time. This means that the proba- 
bility of decay for any one atom is independent 
of its previous history. An atom which has 
survived for a long time has no more chance for 
decay than one formed recently from its mother 
elements. 
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The most convincing proof that the rate of 
decay (and, therefore, also the heat production 
by radioactive elements) has not changed in 
geological epochs can be derived from the study 
of radioactive haloes in minerals. These sharp 
discolorations occurring in concentric spherical 
zones around a radioactive nucleus in certain 
minerals are due to the a-particles emitted 
during long epochs of time. The diameters of 
concentric halo rings correspond very accurately 
to the respective ranges of the a-rays from the 
radioactive nucleus in the mineral and it can be 
established that these ranges, and, therefore, 
the radioactive decay of the imbedded thorium 
and uranium products, do not change within a 
billion years. 

From all these findings we must conclude that 
radioactive elements are giving off heat at the 
same rate whether they are buried deep in the 
mantle of the earth or situated near the surface 
of the earth’s crust and that their heat produc- 
tion was not different at the time when the 
solidification of the lithosphere took place, some 
two billion years ago. 

Chemical and geological considerations indi- 
cate that uranium and thorium may have been 
concentrated mainly in the outermost layers of 
the earth, the lithosphere. R. W. Swanson‘ 
writes: 

“We know that thorium and uranium will rob 
iron of its oxygen, thus they must lie in the 
oxide layer. We know, furthermore, that they 
tend to concentrate in the oxygen-rich rather 
than the oxygen-poor oxides and in fact are 
found most abundantly in the very latest 
magmatic concentrates.” 

“It is to be expected then that the radioactive 
material of the earth will be concentrated in the 
upper layers of the earth.”’ 

The degree in which this concentration of 
radioactive matter actually takes place in the 
direction from the deeper parts of the mantle to 
the surface is unknown. As a matter of fact, we 
would consider it as very improbable that the 
deepest parts of the mantle and even the core of 
the earth are entirely free from radioactive 
substances. This view is supported by results of 
analysis of stone and iron meteorites. The latter, 
which may actually represent the composition 
of the core of our globe, still contain 10-7 gram 
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uranium and 3.10-“ gram radium per gram, at 
least 1/100 of the respective figures for acid 
igneous rocks; stony meteorites contain four 
times as much, per gram, which content is of the 
same order of magnitude as that of the ultrabasic 
deep-seated rocks. 

As far as we can see, the concentration of 
uranium, thorium, and radium with their succes- 
sive products, decreases from acid igneous to 
basic rocks in the ratio of 5 : 1, on the average. 
The total heat generated in these various rocks 
varies accordingly from 5-10~-* to about 1.6-10-® 
calorie per year per gram of rock® where also 
the heat generation of potassium is taken into 
account. 

Arthur Holmes values® are from three to six 
times higher. Evans and Goodman conclude that 
on account of their recent improvements in the 
technique of such determinations ‘‘a substantial 
downward revision is necessary in the present 
conception of the radioactive content of igneous 
rocks. . . . In particular, Holmes’ estimates of 
the rate of generation of radioactive heat 
within the earth may need to be substantially 
decreased.” 

The all-important question for us is now to 
decide whether the higher average values of 
Holmes or the new ones accumulated by Evans 
and Goodman are more representative for the 
computation of the thermal equilibrium of the 
earth. As Slichter® points out in his excellent 
review of the problem of the cooling of the earth, 
in general the figures arrived at by Evans and 
Goodman are only about 1/3 of those which were 
generally accepted within the first quarter of 
our century. 

Evans and his associates certainly had cogent 
reasons for excluding a great number of radium 
and thorium determinations of rocks by previous 
observers as too high. For instance, slight 
contamination of chemicals always tends to give 
too high values of the radium content. There 
are, however, always possibilities of falsification 
of results by influences whose tendency would 
point in the opposite direction. Personally, I do 
not believe that, aside from some inaccuracies of 
the radium standard solutions at that time, 
radium determinations of rocks by Lord Ray- 
leigh (R. J. Strutt), Arthur Holmes, and H. 
Hirschi between 1910 and 1925 were inferior to 
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the most modern ones. Thorium determinations 
naturally are much less reliable on account of 
the lack of equilibrium or of lack of definition of 
the thorium standard solutions used by most 
authors. Direct determinations of uranium are 
probably the least reliable of all and are in most 
cases impossible. The introduction of alpha-ray 
counting methods by Evans and his collabora- 
tors, marks a definite step forward in this respect. 
Nevertheless, I should like to point out a serious 
difficulty which seems to have been overlooked 
so far and which arises when one accepts the 
much lower general mean figures of Evans and 
Goodman for the radioactivity of rocks. 

It is well known that gamma-rays from the 
radioactive substances (in the soil) produce a 
very marked ionization in air or any other gas. 
If-we place an ionization vessel connected to a 
sensitive electrometer, first on the soil or directly 
over a rock, and then over water, we find that 
the observed ionization (expressed in pairs of 
ions produced per cc per sec. of standard air) 
over land is from 2 to 4J greater than over 
water (J=ion pairs per cc per sec.). This ‘‘earth 
radiation” is entirely due to the action of 
gamma-rays of the uranium-radium and the 
thorium family and also of potassium in the soil 
or rock. 

The computation of this effect is very simple 
since we know the absorption coefficients of 
gamma-rays in matter very well. For the 
gamma-rays from radium (RaC), the ‘‘mass 
absorption coefficient” yu/p can be taken as 
4.5X10-* cm?/g. Thus the integrated action of 
gamma-rays from the soil at a point on the 
surface can be taken as:’ 


2rK 
Cana ™=— : 
K/p 





. Mra, 


where gra denotes the ionization (in ion pairs per 
cc per sec.) produced by gamma-rays from 
radium (RaC); K is called Eve’s constant and 
denotes the ionization produced by gamma-rays 
from radium at a point 1 cm from the source. 
K was determined carefully in zine or iron 
vessels of 3-mm thickness as 5.5X10° by the 
author.* Mp, denotes the amount of radium per 
gram of rock or soil. 
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Similar formulas can be used for the computa- 
tion of the gamma-ray ionization of the thorium 
and potassium in the soil. The thorium compo- 
nent of the gamma-rays from most soil and rock 
samples may be taken as almost equal to the 
radium component.’ * The potassium component 
has been overestimated in recent years and a 
revision of the figures is necessary. (A detailed 
account of this analysis will be published 
in Terrestrial Magnetism and Atmospheric Elec- 
tricity, Volume 48.) 

If we adopt Evans and Goodman’s figures for 
the radium content of the four representative 
classes of rocks (acidic igneous, intermediary, 
basic, and sedimentary), we obtain the results 
given in Table I. 


TABLE I. 


Gamma-ray 


Radium content ionization 
Rock type (after Evans and Goodman) Ra 
1. acidic igneous 1.37X10-" g Ra/g 1.05 I 
2. intermediary 0.51 10-" 0.39 I 
3. basic 0.38 X 10-" 0.29 I 
4. sedimentary 0.57 X 107” 0.44 I 


For thorium the respective figures gr, can be 
taken as almost equal, but certainly not higher 
than for radium. 

For potassium (using data cited by Evans 
and Goodman) we obtain Table IT. 


TABLE II. 


Potassium 
content 


Gamma-ray 
ionization 


Values rock type 2K'g qK 
1. acidic igneous 2.8X 107? 0.34 I 
2. intermediary 2.0 10™? 0.25 I 
3. basic 1.4107 0.17 I 
4. sedimentary 2.0 10-2 0.25 I 


Summing up the gamma-ray ionization of all 
three components gra, Grn, GK We obtain the 
y-ray effects given in Table III. 


TaBLe III. 
1. over acidic igneous rocks 2.4 I (qrat+qrnt+ qk) 
2. over intermediary rocks 1.0 J 
3. over basic rocks 0.75 I 
4. over sedimentary rocks is 2 
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The figures given in Table III are maximum 
values, giving the expected ionization due to 
gamma-rays from the radioactive substances in 
the rock. It is obvious that over soil, the effects 
will be considerably smaller on account of the 
admixture of humus and other organic inactive 
decay products. 

Summing up, we may say that if Evans and 
Goodman’s estimates of the average radium, 
thorium, and potassium content of the main 
rock types are correct, we could not expect that 
the difference between the ionization in a closed 
metallic vessel placed over land and over water 
should exceed 1 to 2 J. Actually, much greater 
differences have been observed. Even over sedi- 
mentary ground (gravel), the author actually 
found a difference of 2.77 and many other 
authors report differences of 3-4 J between land 
and water, under similar conditions. 

This discrepancy is serious. I think that direct 
experiments on the gamma-ray effect of certain 
rock or soil types listed in Evans and Goodman's 
report are very necessary and I may add that 
experiments of this type are planned by the 
Department of Terrestrial Magnetism of the 
Carnegie Institution of Washington in collabo- 
ration with Fordham University, upon my 
suggestion.!° 

If powdered rock of a known content (M) of 
radium, thorium, or potassium is placed in a 
concentric shell of known thickness around an 
ionization chamber of spherical shape, the 
gamma-ray effect in the vessel can be computed 
from the formula," 


4nKM 
q= (1 7 e*4) ’ 
u/p 


where p is the density and yu/p the mass absorp- 
tion coefficient of the respective gamma-radiation 
in the layer of. powdered rock (thickness d). 
Thus a direct comparison of results of other 
methods (emanation method) on samples of 
sufficient quantities is possible. Calibration with 
material of standard radium or thorium content 
is relatively easy. Correction for loss of emana- 
tion is small if the granules of rock chosen are 
coarse. 

It is obvious that the above-mentioned dis- 
crepancy between the calculated and observed 
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gamma-ray effects of different soil materials 
disappears if we adopt the higher values of 
radium and thorium contents as given by 
Holmes and by Jeffreys. 

I personally feel inclined to have more 
confidence in these older average values unless 
the discrepancy mentioned above is fully clarified 
or removed in the course of the experiments 
planned. 

Let us now envisage the thermal consequences 
of a change in the adopted figures for the radio- 
activity of rocks. L. B. Slichter® has worked out 
the theory of the cooling of the earth for a great 
number of assumptions in a very admirable and 
complete manner and it is now very easy to 
foresee the final consequences of any change in 
the assumptions with regard to the thermal 
conductivities or the amount of heat developed 
by the radioactive substances embedded in the 
rocks. 

From the radioactive and chemical standpoint, 
the assumption of a layered earth with a rela- 
tively high radioactive heat generation in the 
crust, a lesser one in the mantle, and a much 
smaller one in the core, seems to be justifiable 
and this also suits best the seismological evidence 
of the existence of discontinuities or boundaries 
between these layers. 

In my opinion a decrease of radioactivity with 
depth according to an exponential law, could 
not be upheld as probable for the whole distance 
from the crust to the center of the globe. Never- 
theless, a certain decrease in depth of the 
radioactive matter within each of the three main 
layers may exist. 

Zero radioactivity in the core is extremely 
unlikely since some uranium and thorium must 
have been retained among the heavy elements 
accumulated in the core if the whole earth was 
liquid in the initial state. 

The heat flow derived from the observed 
geothermal gradient varies in different locations 
from 0.6 to 2.5-10-* cal./cm? sec. with an 
average of 1.2-10-*. Thermal conductivities of 
the granitic layer are well known, but for 
greater depth the possible assumptions for the 
thermal conductivity vary very widely. Slichter 
has treated mathematically eleven cases with 
different assumptions as to the radioactivity and 
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thermal conductivity and at different depths 
and arrives at very divergent equilibrium tem- 
peratures ranging from several hundred to many 
thousand degrees centigrade in the mantle and 
in the core of the earth. 

If we exclude some of these cases as most 
improbable, the range of equilibrium tempera- 
tures is considerably narrower. A layered earth 
with radioactive matter decreasing in concen- 
tration with depth from layer to layer but still 
giving in the core about 4 percent of the heat 
production of the granitic layers, seems to be 
most likely (Slichter’s case 3 and 4). Using the 
higher figures of radioactive heat according to 
Holmes and Jeffreys, Slichter computes equi- 
librium temperatures of 450 to 2100°C in the 
core and 1230 to 1360°C in 200-570-km depth. 
This would definitely indicate that the core is 
(at the pressures prevailing there) solid and 
naturally the equilibrium temperatures would 
be still somewhat lowered if Slichter had chosen 
the smaller average figures of radium and thor- 
ium content reported by Evans and Goodman. 

The final or so-called equilibrium temperatures 
for the core and the deeper parts of the mantle 
are still very uncertain; they are subject to very 
large changes if one selects slightly different 
values of the thermal conductivities at great 
depth. A continuous rise of temperature of the 
core on account of the however small radio- 
active heat generated there, is quite possible, 
on account of the poor heat conduction through 
the mantle. As Slichter pointed out, a rise of 
400° in two billion years might be expected if 
the core contains not more radioactive matter 
than the known iron meteorites. Such a process 
would not influence the surface heat flow 
perceptibly. 


SUMMARY 


The heat production of the radioactive sub- 
stances embedded in the rocks is very important 
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for the treatment of the problem of the cooling 
of the earth. Recent average results of Evans, 
Goodman, and others are much lower than 
previously accepted figures. 

It is pointed out that if one computes the 
ionization produced by the gamma-rays from 
radioactive substances in rocks and soil from 
Evans and Goodman's figures, the expected 
ionization turns out considerably smaller than 
the one actually observed by placing ionization 
vessels over land and water. An experimental 
study of this discrepancy is planned. 

It is shown that plausible assumptions as to a 
layered depth distribution of radioactive matter 
lead to temperatures in the core of only 2000°C, 
even if the older and higher data for the radio- 
active heat production in the rocks are accepted. 
A gradual heating of the core itself is to be 
expected (Slichter) if the amount of radioactive 
matter in the core is not smaller than that in 
iron meteorites. 
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The Interior of the Earth 
Viewed in Relation to Earthquake Causes 


III. What Seismology Can Tell Us About the Structure 
of the Earth’s Interior 


By ERNEST A. HODGSON 


Dominion Observatory, Ottawa, Canada 


N accepting a commission to present a paper 

on a prescribed subject, a speaker tacitly 
reserves the right to make his own interpretation 
of the title assigned. Let us examine the limita- 
tions which one might define when considering, 
from the viewpoint of seismology, the symposium 
subject: ‘“The Interior of the Earth Viewed in 
Its Relation to Earthquake Causes.” 

The first limitation which seems justified is 
that the discussion should be confined to that 
part of the earth’s interior where conditions or 
potentials exist which do or which may con- 
tribute to the cause of earthquakes. As most 
earthquakes originate at depths ranging from, 
say, a kilometer to about 700 km, one might at 
first thought decide to consider only the spherical 
shell within the earth which is bounded by these 
two horizons. 

That there is ample scope for discussion even 
with this stricture may be taken for granted. 
It would be desirable, however, to extend the 
depth range to include, for secondary considera- 
tion, those horizons in which conditions exist 
which contribute to the data on which are based 
the accepted the 
earthquakes. 


theories as to cause of 

If this be granted, the upper limit must be 
raised to the actual surface of the earth and the 
lower set at its very center. For the mountains 
and the ocean basins, the geosynclines and the 
tectonic structures of the outer crustal layers 
are certainly surface indications which contribute 
to our deductions as to the cause of shallow-focus 
earthquakes. And the discontinuities, especially 
the great discontinuity at the core, together 
with the distribution of elastic properties and 
density give similar evidence on which are based 
the accepted conclusions as to the cause, and 
even the very existence, of deep foci. 
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It may be assumed, then, that the primary 
concern in this presentation is that range of 
depths in which earthquakes originate but that 
it is proper to admit the discussion of so much 
of the rest of the earth’s interior as may con- 
tribute to an understanding of the cause of 
earthquakes. Further, since the viewpoint is that 
of seismology, the paper is limited to an exami- 
nation of those data which have been obtained 
directly from studies of earthquakes and of their 
instrumental records. 

Dutton! expressed what was the accepted view 
of seismologists of his day when he set the 
maximum depth of earthquake foci at twenty 
miles. When, in 1921, Walker* deduced probable 
depths as great as 800 miles, his conclusion did 
not find general acceptance. It was supported, 
however, by a series of papers by Turner,* of 
which the first appeared in 1922. Wadati’s' first 
paper on the distinctions indicating the existence 
of shallow and deep-focus earthquakes appeared 
in 1928. 

Since then, the evidence has grown from the 
studies of many, so that today seismologists 
generally accept the conclusion that earthquake 
foci may appear at depths as great as 700 km. 
One of the earlier papers, which served to con- 
vince even the most reluctant, was the doctorate 
thesis of Father Stechschulte,® published in 1932. 
A paper by Scrase,® based on studies carried out 
about the same time but published a few months 
earlier, set forth in some detail the characteristics 
of deep-focus earthquakes as shown by seismo- 
graph records. The graphs, published by Father 
Brunner’ in 1935, brought the means for ready 
evaluation of focal depths from seismograms 
into the hands of station operators, who were 
thereby able to prove readily for themselves, 
through the evidence shown by their own current 
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records, the inescapable conclusion that earth- _ 


quakes do originate in some cases at very 
considerable depths, or ‘‘finite depth’ as Walker 
put it two decades since. Father Macelwane® 
gave the name “plutonic earthquakes” to this 
newly accepted class of foci. Within the past 
few years, a series of papers by Gutenberg and 
Richter® has brought together an overwhelming 
mass of evidence concerning them, so that today 
any hypothesis of earth structure in the zone of 
earthquake foci must take cognizance of these 
modern data. 

Accepting for the moment the statement that 
earthquakes originate at foci which may be 
classified as ‘“‘shallow’’ and “‘deep,”’ attention 
may be directed to the theories advanced to 
explain the cause of earthquakes. We shall 
consider initially only shallow-focus shocks, the 
so-called ‘‘normal” earthquakes, and may profit- 
ably pass over the early philosophical theories as 
to their cause, which range from the supernatural 
to the pseudoscientific. One of the first inescap- 
able deductions from the data obtained through 
the world-wide network of seismograph stations 
set up through the efforts of Milne, toward the 
end of the nineteenth century, was that they 
could not be generally due to volcanic action. 
The natural division into volcanic and tectonic, 
or structure building earthquakes followed. It is 
with the cause of tectonic earthquakes that we 
here concern ourselves. 

As observed by Leet,'® there is ample evidence 
to show that whatever be the origin of the forces 
at work to change the structure of the earth “‘it 
may be accepted as axiomatic that the cause of 
earthquakes is intimately related to the universal 
mobility of the earth’s surface."’ Some sugges- 
tions which have been offered as to the mecha- 
nisms by which these forces might be built up are: 
the shrinkage of the earth due to cooling, the 
variation with latitude and depth of the hori- 
zontal component of centrifugal force (‘‘Pol- 
‘fluchtkraft’’), the transfer of sediments, and the 
removal of great ice sheets. Certainly the cause 
of shallow-focus earthquakes is closely related to 
orogenic or epeirogenic forces. Critical exami- 
nation of the above mechanisms by many 
specialists shows that each seems insufficient, by 
itself, to account for the evidences found by 
geologists on the surface of the earth. Each 
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implies a rigid crustal layer of considerable 
strength overlying one, more or less deep, of 
weak, but dense and viscous, material. This 
concept is supported by the measurements of 
gravity and geodesy. We confine our attention 
to the evidence from earthquakes and their 
records. 

In 1906 a great earthquake displaced the San 
Andreas fault in California throughout a surface 
span of more than 200 miles. Some horizontal 
displacements of adjacent edges of this fault were 
more than twenty feet. After studying the 
evidence obtained, Reid" proposed the “‘elastic- 
rebound” theory to account for the data. 
According to this theory the forces, in this case 
largely horizontal, had been built up by some 
one or more of the mechanisms mentioned or by 
others not yet suspected. These forces tending 
to cause the shift were opposed by friction, thus 
storing up in the distorted edges of the fault the 
energy which was released at the time of the 
earthquake. This explanation is generally ac- 
cepted in the case of normal earthquakes. 

That some earthquakes result from forces with 
a strong vertical component is evidenced by 
many field studies. The vertical displacement in 
the Hawke Bay earthquake in New Zealand in 
1931 was, in places, as great as eight feet, while 
in the Alaskan earthquake of 1899 the sudden 
vertical uplift was nearly fifty feet. It may be 
granted from field studies of shallow-focus earth- 
quakes alone that both horizontal and vertical 
movements occur with many examples of a 
combination of these, resulting in oblique shifts 
at the exposed fault scarp. In each of these, 
neglecting the unsettled question of the mecha- 
nism which caused the strain to develop, it may 
be taken as a fact that the earthquake was a 
manifestation of a release of strain as set forth 
in Reid’s theory of elastic rebound. 

Other papers in this symposium deal with the 
subject of isostasy. The only observation to be 
made here is that the depth of isostatic compen- 
sation is not likely to be a discontinuity as 
evidenced by the travel times from earthquakes. 
Daly" states that the elastic properties and the 
density of a rock changing from crystalline to 
vitreous conditions at depth vary in about the 
same ratio, so that a sudden transition in seismic 
wave velocity is not to be expected; and, 
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according to Gutenberg," it cannot be observed. 
We are not here concerned directly with wave 
velocities except insofar as they reveal condi- 
tions which have some bearing on the cause of 
earthquakes. The Mohorovitié discontinuity, at 
a depth of 30-40 km is, according to general 
seismological opinion, not to be identified with 
the depth of isostatic compensation at a depth 
variously estimated at from 10-120 km but 
certainly well within the asthenosphere. Not 
more than two decades past it was generally 
supposed that earthquake foci could not occur 
in the asthenosphere. 

The seismological evidence that foci do occur 
at such depths is overwhelmingly conclusive. 
Walker deduced the fact from computed angles 
of emergence, based on data furnished by 
Galitzin. The attention of Turner was first 
directed to it by the too-early arrival times at 
the antipodes of seismic waves from certain foci 
which he defined as “‘deep.’’ Two phases due to 
longitudinal waves reflected at the surface should 
occur if the focus is appreciably deep, as was 
pointed out by Walker? as early as 1921. These 
two phases have been identified in many studies 
of records from plutonic earthquakes. There 
should be a marked reversal in curvature in the 
time-distance curves for deep foci. A number of 
earthquakes have been intensively studied from 
the world records, as was that of the Japanese 
earthquake of March 29, 1928, by Stechschulte.® 
The reversal of the time-distance curve has thus 
been firmly established. 

Certain characteristics identify a record at 
once as having been due to waves from a deep 
focus. The phases are more numerous and are 
sharply defined. There is a markedly smaller 
surface phase. The interval between the arrival 
of the direct longitudinal and transverse waves 
is abnormally long for short epicentral distances. 

The characteristics of earthquakes from deep 
foci are no less pronounced than those of their 
seismograms. Although the latter prove the 
shock to have been very severe, the surface 
disturbance, even at the epicenter, is relatively 
small. Wider areas feel the tremors more than is 
the case for shallow-focus shocks. The direction 
of the motion at or near the epicenter is notice- 
ably vertical. 
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The above paragraphs do not, by any means, 
complete the list of evidence proving the occur- 
rence of foci at considerable depths. It is not the 
province of the present paper to detail these nor 
to outline the various means for determining the 
depth. The above sources of evidence have been 
given as a brief indication to those approaching 
this subject (of this symposium) from other 
fields of research that the fact of deep foci has 
been established beyond doubt and that the 
determinations of epicenters and focal depths for 
a large number of these have been carried out 
with at least fair accuracy. Some, studied in 
great detail, serve to check the rougher estimates 
made for others. A great mass of data has been 
gathered in the two papers by Gutenberg and 
Richter entitled, ‘‘Materials for the study of 
deep-focus earthquakes.’’® Similar lists were 
published earlier by Turner,* Wadati,‘ Berlage," 
Stoneley,'® Sharpe,'*® Lynch,” and others. 

We are here concerned with the geographical 
distribution of the epicenters of earthquakes and 
with their focal depths, for these are inescapably 
related on the one hand to the cause of earth- 
quakes and, on the other, contribute to our 
knowledge of the interior of the earth. But before 
discussing this ultimate phase of the subject we 
must outline one further pertinent phase of 
seismological evidence. 

In some earthquakes the direction of the 
initial movement of the ground at a distant 
station, due to the arrival of the first longitudinal 
wave, is down or toward the epicenter, in others 
it is up or away from the epicenter. Movements 
of the first type are called kataseismic, those of 
the second type anaseismic. It has been found 
that earthquakes show anaseismic motion at 
some stations and kataseismic motion at others. 
Further study has revealed other facts about 
this distribution which are exceedingly interest- 
ing but which do not contribute to the present 
subject. One feature does do so, however. It has 
been found that initial motions at any given 
station for all foci in a given epicentral region are 
almost always either kataseismic or anaseismic. 
The implication of this observation is far 
reaching. 

It was early suggested that deep-focus earth- 
quakes might be due to explosions at depths so 
great that the lack of strength in the rocks 
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would, it was argued, preclude the possibility of 
slippage along fault planes. If that were the 
case, then all such deep-focus earthquakes should 
always record as anaseisms. As they do not, the 
motion being kataseismic at stations in some 
azimuths from the epicenter and anaseismic in 
the others, we are almost forced to conclude that 
the mechanism of plutonic earthquakes is the 
same as for those originating at lesser depths, 
i.e., slippage on a fault plane. The ‘strong 
transverse waves generated from deep foci are 
further evidence supporting this assumption. 

We come now to the distribution of earthquake 
epicenters and their foci. Before any distinction 
was made between the two classes of foci, many 
earthquake maps were constructed. These showed 
with increasing emphasis, as the more accurate 
seismological data accumulated, that, in general, 
the regions of greatest seismic activity were 
closely associated with mountain building. A 
belt of high seismicity circles the Pacific. Another 
crosses through the Dutch East Indies across 
the north of India to the Mediterranean. A 
third, of minor importance, begins near Iceland 
and traverses the mid-Atlantic. Also, as time 
went on, a further truth emerged with consider- 
able emphasis—severe earthquakes sometimes 
occur quite outside these belts. 

We now speak of earthquake foci at depths up 
to about 60 km as “normal,” those in the 
approximate range 60-250 km as “intermediate,” 
and those at 250 km and greater as “deep.” 
While normal and intermediate earthquakes may 
occur outside the seismic belts, deep foci occur 
only in those belts and only in five general 
regions thereof, namely: Japan, the Dutch East 
Indies, New Guinea, the New Hebrides, and in 
the western part of South America, respectively. 
We may examine in some detail the Japanese 
group. 

In his earlier papers, Wadati concluded that 
it was possible to draw lines of equal focal depth 
through the epicenters of earthquakes plotted 
on a map of the Japanese islands. These lines, 
which might be called isobathyseismals, were, 
as he drew them, wide arcs roughly parallel to 
the western edge of the Kurile-Japanese deep. 
Later, Wadati abandoned the isobathyseismals 
but defined two deep-focus zones. One of these 
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lies approximately northwest-southeast or trans- 
verse to the main Japanese island, Honsyu, a 
little south of Tokyo, and extends up into the 
Sea of Japan and down into the Pacific. A 
second zone, of lesser importance, runs almost 
at right angles to the first from its upper end and 
bears roughly northeast across the northern tip 
of the second largest island, Hokkaido. 

The plotting, by Gutenberg and Richter,® of 
the locations, carefully computed from all the 
reliable data available, confirm these two zones 
defined by Wadati but reveal a striking fact: 
that the deep foci all lie on the edge of the zones 
farthest from the Pacific. This same arrangement 
is found for the Dutch East Indies region and 
for South America where the deep foci are found 
to the east of the Andes. The shallow foci lie 
nearest the ocean deeps and the intermediate 
ones between the other two. So far, the South 
American data reveal no foci between 290 and 
600 km. 

It must be noted that data on depth of focus 
are not available for more than thirty years past 
and that only during the last decade have they 
been of present-day accuracy. It is, therefore, 
not very safe to generalize. However, for what 
it is worth, we may sum up by saying: that the 
fact of deep foci is well established; that various 
criteria for determining depth are known and 
yield concordant values which are of high 
percentage accuracy; that the deep foci of which 
we so far have record are confined to the five 
regions mentioned; that they are associated with 
intermediate and shallow foci in a distribution 


- which roughly shows the deep foci farthest from 


the oceanic deeps and the shallow foci nearest; 
that the deep foci show anaseismic and kata- 
seismic properties which are generally the same 
as for intermediate and shallow foci in the same 
region; and that deep foci show strong transverse 
waves on their seismograms. 

The interpretation placed on the above is that, 
in the five regions mentioned, there are deep 
oblique, thrust faults which we might call geo- 
fault-planes and that all the foci result from 
slippage on these great planes. 

Accepting this interpretation, we may examine 
a little more closely the phenomena of deep foci. 
It has been found that these “‘repeat’’ as do 
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shallow foci. That is to say, shocks of almost 
equal intensity are recorded from the same focus 
at intervals of two years and less. If the interpre- 
tation of the paragraph next above be true, then 
forces must be at work building up strain at 
depth at a remarkably high rate. Moreover, the 
fact of repeat shocks is a mortal blow to the 
theory of an explosive cause for these deep foci. 

It is also found that aftershocks occur in 
connection with deep foci as they do for those 
more shallow, again indicating a common 
mechanism. They are probably quite as numer- 
ous, but the smaller ones fail to record at the 
great distances due to the depth involved, just 
as small aftershocks at normal foci record only 
at near stations. 

In certain regions the horizons at which the 
deep foci occur seem to be restricted. For 
example, in South America, many foci occur at 
depths from 600 to 660 km with epicenters east 
of the Andes, while intermediate shocks at 
depths of 100 to 290 km have epicenters in the 
Andean zone with strong normal shocks in the 
oceanic deeps and the coastal plane. In other 
regions the distribution in depth is almost 


continuous, as in Japan and in the Dutch East 
Indies. 

A strong deep focus earthquake occurs on the 
average about once a year. There seems to be no 
particular relation between their intensity and 
their depth, some of the strongest being at 
nearly 700 km. So far, again we repeat, so far, 
no depth appreciably greater than 700 km has 
been found. 

This seems to present the case from the 
seismological viewpoint. It is perhaps not out 
of place, in conclusion, to point out that the 
laboratory experiments of Griggs, Bridgman, 
and others show that plastic flow of rocks at 
high confining pressures and temperatures ex- 
hibits a stop-and-go phenomenon which lends 
strength to the deduction that, even in the 
asthenosphere, the viscosity is so high that, if 
the forces build up with sufficient rapidity, 
slippage along fault planes could occur as in the 
lithosphere. The repeat foci show that the 
building up of these forces is remarkedly rapid. 
The question of the mechanism by which they 
so accumulate lies beyond the scope of the 
present discussion. 
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HAT is it that makes the “pure” physicist go when 

he is on the trail of some new idea in his laboratory ? 
The answer is, of course, complex, but I believe that 
through all the multifariousness runs one simple guiding 
thread, the craving for understanding. To the extent that 
the guiding motif of an enterprise is the craving for under- 
standing, to that extent the enterprise may be said to be 
purely scientific, as distinguished from technological, or 
utilitarian, or artistic, or political, or what not. The craving 
for understanding reaches its greatest poignancy only in 
a few cases, but all of us who are engaged in pure research 
have it to a certain extent, and it is the vital part of what 
makes us go. It is not a matter to be argued about, as to 
whether such a craving has economic or other justification; 
it is only to be accepted as a fundamental fact about human 
beings that some of them have developed to a high degree 
the passion for understanding and a delight in the corre- 
sponding activities, just as others have a strongly developed 
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sense of beauty or of conduct. If society is ever going to 
become anything more than a vicious merry-go-round of 
circular activity, if ever there are ends in themselves or 
goods in themselves, then surely the gratification of the 
craving for understanding is one of them. 

To those who have a passion for understanding, society 
will not be a satisfactory place unless it affords opportunity 
for the acquiring of understanding, so that to the extent 
to which the function of society is to make life satisfactory 
for its members, and it seems to me that this is pretty 
nearly the whole function of society, one of its responsi- 
bilities is the making and providing of adequate scientific 
opportunity. Society is the servant of science even more 
and in a more fundamental sense than is science the servant 
of society. Any control which society exerts over science 
and invention must be subject to this condition.—P. W. 
BRIDGMAN, in “Science, and Its Changing Social Environ- 
ment,”’ Science, January 12, 1943. 


JOURNAL OF APPLIED PHYSICS 

















The Interior of the Earth 
Viewed in Relation to Earthquake Causes 


IV. Measurements of Terrestrial Magnetism Tell Us Little 
About Earthquakes but May Tell Us Much 
About the Interior of the Earth 
By WILLIAM A. LYNCH 


Fordham University, New York, New York 


EARTHQUAKES AND MAGNETIC EFFECTS 


N the old treatises on earthquakes, the authors 

were careful to describe at great length the 
other physical phenomena that seemed to be an 
accompaniment of the seismic disturbance; the 
most common of these were lightning, auroral 
displays, and local magnetic effects. For example, 
Ponton in Earthquakes and Volcanoes, published 
in 1871, mentions repeatedly that ‘“‘magnetic 
needles are strongly affected and permanent 
magnets drop their keepers ;”’ he gives as a theory 
of the cause of some “mild earthquakes,” the 
discharge of electricity from the air to a cavity 
within the earth; as a proof he cites the circum- 
stance that ‘in such cases, permanent magnets 
lose their magnetism for a short time, and drop 
their keepers before the shock is felt.”’ 

Knott in his text The Physics of Earthquake 
Phenomena, published in 1908, summarized the 
comprehensive investigations by Milne on con- 
nections between earthquakes and magnetism by 
saying “‘the only certain conclusion that can be 
drawn is that when the seismic disturbance is 
large enough, the magnetographs act to some 
extent like seismometers. In other words, the 
effects were dynamical, not magnetic. In the 
India earthquake of April 4, 1905, quite distinct 
disturbances were observed on the magneto- 
graphs at Bombay; but an inspection of these 
quite bears out the statement of the director, 
Mr. Moos, that they are mechanical effects.” 

Chapman and Bartels in their two massive 
volumes on Geomagnetism mention earthquakes 
only once, to this effect: ‘“‘Earthquakes cause 
mechanical oscillations of the magnetic needles, 
but no certain magnetic effect has so far been 
found.” 

McNish in 1937 at a symposium in Washing- 
ton on the ‘Interior of the Earth’’ dismissed the 
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possibility of a connection existing between 
terrestrial magnetism and earthquakes with the 
statement: ‘Efforts to discover magnetic changes 
attributable to earthquakes have met with no 
success.” In Earthquakes by Captain Heck, mag- 
netician and seismologist, the word magnetism is 
not mentioned, except on the title page, where 
the official position of the author is given. 


ENERGY RELATIONS 


Thus it would appear that my task is an easy 
one in this discussion; I can with complete 
justification state that magnetism can be ignored 
as a cause of earthquakes and close this part of 
the paper. But the question is worthy of more 
consideration ; it seems worth while to investigate 
the possibilities that exist and dispose of them 
one by one. First consider the matter from the 
viewpoint of the energy associated with the 
magnetic substance and field of the earth. Ac- 
cording to Bauer, the magnetic properties of the 
earth are due largely to sources in the interior, 
only about 6 percent of the total magnetization 
being ascribed to external causes; for the pur- 
poses of this part of the discussion, I shall ignore 
the external causes and confine my attention 
solely to the origin of the main field. Various 
models can be used to describe the main field; 
we may assume a single magnetic dipole at the 
earth’s center, the moment of the dipole being 
M = H,a* =8.1 X10” gauss cm’. (4/7) =0.315 gauss 
=field at magnetic equator, Chapman and 
Bartels, p. 645.) Alternatively we may assume 
that the earth is a uniformly magnetized sphere 
of intensity of magnetization J=0.08 gauss (this 
is many orders of magnitude greater than is 
observed in ordinary rocks) ; the external fields 
are the same in the two models. In order to 
visualize the second model, we may assume with 
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Chapman and Bartels that we have parallel 
cobalt steel bars imbedded in a non-magnetic 
earth, one bar to every cubic meter of the earth, 
each bar containing 1.5 kg of cobalt steel and 
occupying 200 cc and having /=400 gauss; the 
resulting field will be approximately the same as 
that of the earth. Or we may take McNish’s 
model; he assumes that a highly magnetized 
embroidery needle is imbedded in every 5 cc of 
the earth’s substance. Associated with any vol- 
ume distribution of magnetized matter is a 
definite energy, which amounts to }.M?/a* for a 
uniformly magnetized sphere, where M is the 
equivalent dipole moment (8.1 X 10%) and a is the 
radius (6.37 10° cm for the earth); the energy 
for the whole earth is 1.310” ergs. The energy 
of the Charleston earthquake was about 1.4 x 10 
ergs. Captain Heck estimates the potential 
energy stored before release in the San Francisco 
earthquake as 1.7 X 10" ergs; a calculation of the 
seismic energy released during the slight earth- 
quake in 1938 in the vicinity of Freehold, New 
Jersey shows that it amounted to about 10" ergs. 
Thus we see that the total energy of the earth 
bound up in its magnetization is only that of a 
severe earthquake; in the latter case, the energy 
is drawn from a limited region and the energy 
density must be many orders of magnitude 
greater. No process is known by which the energy 
of magnetization of the earth can be transformed 
into mechanical energy, particularly with the 
suddenness characteristic of earthquake phe- 
nomena. In the Barkhausen effect in iron, we do 
have elementary molecular magnets changing 
their orientation in groups of varying size as the 
external magnetizing field is altered, the process 
being discontinuous; each change in position is 
accompanied by an alteration in the surrounding 
magnetic field and is made evident by an induced 
current in an adjacent coil connected to an 
amplifier, and ultimately by a sound in a loud- 
speaker. Here then is a sudden transformation 
from magnetic energy to mechanical energy, but 
it is difficult to see any other similarity between 
this “‘ironquake”’ and an earthquake. A possible 
mechanism is available by which magnetic energy 
can be replaced if it be transformed into mechani- 
cal energy; Hulburt states that energy flows to 
the earth from the sun during a magnetic storm, 
to the extent of about 3X10” ergs per storm; 
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thus many storms would be required to make up 
for the loss experienced during one earthquake. 
Actually the excess magnetic storm energy is 
dissipated as heat in the earth due to 77 losses 
of induced currents. These various energy values 
are dwarfed to insignificance when compared to 
other aspects of the energy of the earth: The 
energy of rotation alone is about 10** ergs, and 
the energy intercepted from the sun as sunlight 
is 2X 10* ergs every second. 


EARTH ELECTRIC CURRENTS 


A different way of visualizing the earth’s field 
is in terms of surface electric currents circulating 
from east to west, the current density in any 
region depending on the sine of the magnetic co- 
latitude, and the total current being of the order 
of 1000 million amperes; the magnetic field of 
this current would be approximately the same as 
that of the main field of the earth. The magnetic 
energy associated with current flowing in a 
circuit is given by $7? where L is the effective 
inductance of the shell of the earth used as the 
conductor and the total energy is the same as 
that of the uniformly magnetized spherical earth, 
1.310% ergs. The advantage of this picture is 
that we can calculate the time constant of the 
inductive circuit and hence obtain an estimate of 
the time required for the discharge of the circuit 
and the liberation of the energy; but the dis- 
advantages are very much greater for we cannot 
conceive of the shell acting like an ordinary 
circuit which can be closed and opened at will, 
liberating energy at specified points and in 
specified amounts. Furthermore, there is no 
experimental evidence toindicate that surface cur- 
rents of the necessary density exist ; the currents 
measured by the earth electricians are expressed 
in terms of the drop in potential produced per 
km; in spite of the high resistivity of the surface 
layers, the drop is only a few millivolts per km. 
I feel that no useful information can be derived 
by pursuing this aspect of the matter further. 

It seems clear, therefore, that terrestrial mag- 
netism cannot be considered in any sense as a 
primary cause of earthquakes—the necessary 
energy is not present nor is any agency known 
by which the transformation from magnetic to 
mechanical energy can be accomplished even if 
the energy were available. 
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MAGNETIC ACTIONS AS TRIGGER FORCES 


But we may obtain information from terrestrial 
magnetism about trigger forces and the interior 
of the earth. A magnetized body is in a state of 
mechanical strain; the term magnetostriction is 
used to describe all the phenomena that come 
under this classification; and magnetostrictive 
effects have to do with dimensional changes in- 
volving length, volume, and twist. The Joule 
effect is perhaps the simplest of these and is 
made evident by an increase in length of an iron 
rod when the rod is magnetized ; it was discovered 
by Joule in 1847. For iron at low fields, the Joule 
coefficient amounts to about 3X10-* cm/cm/ 
gauss; the coefficient of linear expansion of iron 
is 1.2X10-° cm/cm/degree C. It is clear that the 
clastic stresses due to magnetostriction are insig- 
nificant in comparison with those due to tem- 
perature change and it is unnecessary to pursue 
this topic any further. 

A conductor carrying current in an external 
magnetic field is subjected to a mechanical force 
which is proportional to the current and the 
normal component of the field. Earth currents in 
general are very weak and so too is the earth’s 
magnetic field ; hence the forces of interaction are 
exceedingly small. The maximum current that 
flows due to natural causes occurs in thunder- 
storms and Schonland gives the value of 100,000 
amperes for a single powerful stroke of lightning ; 
the force per cm on such a current filament in a 
field of one gauss is 10* dynes, equivalent to 
about 3 ounce. Thus if the stroke were to pene- 
trate deeply into the earth, which it does not, 
the force on the conducting earth due to the 
magnetic effect would be very small compared 
to the forces of expansion produced by the heat- 
ing effect of the current; the energy of a single 
lightning stroke is about 10" ergs, equivalent to 
that of a minor earthquake. 

The mechanical forces arising within the earth 
due to magnetic effects are very small and can 
be ignored even as trigger forces. 


TEMPERATURES IN THE INTERIOR AND 
THE CURIE POINT 


Terrestrial magnetism may give information 
about the interior of the earth. The concept of a 
uniformly magnetized sphere may be all wrong 
because of the effect of increasing temperature on 
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magnetizable bodies. As the temperature is 
raised, ferromagnetic substances lose their mag- 
netization at the inversion temperature or Curie 
point; this occurs at 780°C for iron, 580°C for 
magnetite, and 350°C for nickel. Near the sur- 
face of the earth, the temperature increases as 
we proceed downward at the rate of 30°C/km; 
it is considered certain that the increase continues 
for some distance below the surface. At 25 km 
the temperature to be expected is that of the 
Curie point for iron and at greater depths all 
ferromagnetic substances may lose their mag- 
netic properties. It is possible that the Curie 
point may be raised with increasing pressure and 
that the change does not occur until a very much 
greater depth is reached. Experiments of Adams 
and Green show that the change in Curie point 
with pressure is less than 10°C for a pressure of 
3600 atmospheres, corresponding to a depth of 
13 km under the continents; they state “‘if there 
is any trend at all, it is in the nature of a slight 
decrease in temperature with increasing pres- 
sure.”’ If the magnetization of the earth is con- 
fined to an outer shell 20 km thick, the intensity 
of magnetization within the shell would have to 
be of the order of 8 gauss; and if the thickness is 
taken as 84 km the magnetization would be 2 
gauss. The greatest known magnetic anomalies 
are those of Kursk and Lapland, produced by 
magnetite with J=less than one gauss; it is 
practically impossible, therefore, for such great 
average intensities as 8 or 2 to occur in the crust. 

Goranson at the Geophysical Laboratory has 
recently attained pressures of 200,000 atmos- 
pheres; a test on the change of the Curie point 
of a CdMgFe spinel gave a shift of 0.0005— 
0.002°C per atmosphere with pressures to 10,000 
atmospheres. McNish states ‘“‘To permit ferro- 
magnetism at the temperatures which are as- 
signed to the earth’s interior, the Curie point 
temperature would have to increase about 0.01° 
for a pressure increase of one atmosphere.”’ 
However, he thinks it possible that the estimated 
rate of increase of temperature with depth is too 
great and should be of the order of 10°/km; 
he is supported in this view by the measurements 
of Van Orstrand on deep wells in the United 
States and South Africa. This rate of increase of 
temperature with depth gives a value of about 
2200°C for the earth’s center; a rise of 0.001°/ 
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atmosphere in the Curie point would permit the 
core to be ferromagnetic. Gutenberg arrives at 
a value of about 2000°C for the center of the 
earth; he bases his argument on the variation of 
temperature with the depth at which the transi- 
tion from the crystalline to the viscous state 
occurs. Hence if Goranson’s findings on the shift 
of the Curie point with pressure are verified, it 
may be that the core, postulated by the seis- 
mologist, is a ferromagnetic mixture of nickel 
and iron. 


CONDUCTIVITY OF THE EARTH AND 
DEEP-FOCUS EARTHQUAKES 

One other aspect of terrestrial magnetism may 
lead to important information about the interior 
of the earth. A daily variation in the earth’s 
field occurs, caused largely by ionization of the 
atmosphere by the sun; ionization currents flow 
above the earth and induce currents in the earth ; 
analysis of the daily variation of the field yields 
the components supplied by each set of currents, 
the currents themselves, and finally the con- 
ductivity of the interior of the earth. The results 


show that the conductivity is low until a depth . 


of 200-250 km is reached, where a great increase 
occurs, from about 10~-" c.g.s. for dry earth or 
rock to 3.6X10-". The daily variation does not 
supply information about greater depths because 
of the shielding effect of the induced field on the 
inducing field. But the magnetic field produced 
by a magnetic storm penetrates to a much greater 
depth ; analysis of the resultant field indicates a 
conductivity of about 4X10-" c.g.s. at a depth 
of 700 km; Lahiri and Price estimate the value 
of the conductivity as about 10~" c.g.s. below 
this depth and suggest that a change occurs to 
a greater metallic content at this depth. Their 
“core’’ would be quite different from that of the 
seismologists. It may be remarked that many 
deep-focus earthquakes occur at depths of 200— 
250 km, especially in the Hindu Kush Moun- 
tains and in Chile; the deepest earthquakes re- 
corded have a depth of focus of about 700 km. 


THE EARTH’S CORE 


Fleming has reported very recently that_ 


MeNish has worked out a system of dipoles to 
represent the residual field of the earth ; the main 
field is represented by the magnetic dipole of 
moment 8.110” c.g.s. units as before, but the 
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residual field is represented by 14 radially di- 
rected dipoles midway between the center and 
the surface of the earth, i.e., at the surface of the 
seismologists’ core; the depth was not chosen 
arbitrarily; a greater or less depth would have 
required a much greater number of dipoles for 
the same closeness of representation. Secular 
change consists almost entirely of changes in the 
residual field; magneticians may hope to find in 
these results a solution of the problem of the 
earth’s field and its secular change. 

Among the most recent suggestions of the 
origin of the earth’s field is that of Elsasser; he 
supposes a circulation within the metallic liquid 
core of the earth which gives rise to a tempera- 
ture gradient in the east-west plane; a net 
thermoelectric current flows encircling the axis 
of the earth and setting up the observed mag- 
netic field. If this hypothesis survives, and doubt 
has been cast upon it, it will not affect any con- 
clusions drawn thus far about the cause of earth- 
quakes; but it will strengthen the evidence in 
favor of a liquid core. 

In conclusion, I may summarize the state- 
ments above by saying: (1) terrestrial magnetism 
does not yield information of importance to the 
seismologist on the question of the cause of 
earthquakes; (2) it may supply evidence of im- 
portance, corroborating the presence of discon- 
tinuities within the earth found by the seismolo- 
gist, particularly that of the earth’s central core. 
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The Interior of the Earth 
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VY. Tectonophysics or the Physics of Earth Deformation 


BY JAMES B. MACELWANE, S.J. 


St. Louis University, St. Louis, Missouri 


ECTONOPHYSICS may be defined as the 

physics of earth deformation. As a science 
in its own right it is very young indeed. Tectono- 
physics deals with the physical properties of 
rocks and minerals as, for example, with their 
densities, elasticities, strengths, internal friction, 
viscosities, crystallographic lattices, tendencies 
to harden or to soften under plastic flow, their 
characteristic resistance to rupture and flow, 
polymorphic transitions under changing condi- 
tions of pressure and temperature, with their 
melting and boiling points, their latent heats, 
their thermal conductivities and diffusivities, 
their exothermic and endothermic reactions in 
so far as these are related in any way to the 
deformation of those rocks or minerals. These 
properties may be and have been studied in 
various ways. Experiments of many kinds have 
been conducted on specially prepared samples 
in the laboratory. Metallurgical methods have 
been applied to rock melts.! Scale models of 
various types have been designed and used.’ 
Some of the above properties have been meas- 
ured on a large scale in the field.*? However, all 
these processes can be applied only to those 
parts of the earth that are accessible through 
exposures, mines, drill cores, and wells. The 
methods of tectonophysics, therefore, are limited 
at present to the upper three miles of the earth’s 
crust. By means of high pressure and high 
temperature apparatus the conditions may be 
simulated for a few decades of miles of depth. 
Beyond this practical limit tectonophysics must 
needs work by inference. In spite of these severe 
limitations, however, the new science of tectono- 
physics already has added much to our knowl- 
edge of the behavior of rock masses and of the 
mechanics of their deformation. The older 
experimental results of F. D. Adams,‘ L. H. 
Adams,® and _ others, the 


and newer results 
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obtained by Bridgman,’ Griggs,’ Birch,’ and 
others on the deformation of rocks and minerals 
have taught us the futility of vague speculation 
concerning the causes and mechanism of large 
scale deformations, the boundary conditions of 
which we know little or nothing; on the jointing, 
folding, faulting, intrusion, extrusion, and miner- 
alization, the end result of which alone is 
observed by the geologists; and of extrapolating 
beyond any of the data we possess far into the 
depths of the earth. These experiments have also 
taught us the necessity of determining for each 
of the materials which actually constitute the 
earth’s crust their specific physical and chemical 
behavior when subjected to large compressive 
and shearing forces. It has been found impossible 
to generalize. Each substance behaves in its own 
characteristic fashion. A lubricant like graphite, 
for example, was found by Bridgman to cease 
flow and to become an abrasive when sheared 
sufficiently. Most substances deform by jumps 
and refuse eventually to flow at all. Silica glass 
when subjected to shearing forces will not flow 
but instead tends to become more like a quartz 
crystal. In general Bridgman found that the 
resistance to shearing strain increases with the 
hydrostatic pressure exerted on the material. 
In regard to the length of time of application of 
the deforming forces, each substance seems to 
have its own pattern of behavior. Bridgman, 
for example, found that some materials decrease 
their resistance appreciably as the speed of the 
deformation is decreased, while mica and the 
larger number of substances showed no meas- 
urable variation of force with speed of deforma- 
tion. When a constant force is applied some 
substances as, for example, silver, begin to flow 
rapidly then gradually slow down and eventually 
refuse to flow, while others like mica become 
weaker so that the rate of yield increases until 
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instability is reached and the substance becomes 
a lubricant. It is clear from these experimental 
results that the method of scale models must be 
used with extreme caution. How a given geo- 
logical structure will behave when subjected to 
a specified force system or field of force will 
depend not only on the boundary conditions, on 
the geometry of the structure, and on the 
physical and chemical environment in which it 
is placed, but also on the characteristic behavior 
of each individual constituent under these 
conditions. It is also clear from the experimental 
results that we have no means at present of 
predicting the behavior of a substance beyond 
the application of certain more or less obvious 
considerations such as orientation with respect 
to cleavage planes, type of crystallographic lat- 
tice, molecular structure, and similar properties. 

It is clear then that the theory of elasticity 
will not suffice for the analysis of earthquake 
phenomena. Very much must be borrowed from 
hydrodynamics. But the other extreme will not 
work either. The earth cannot be treated as a 
fluid body. We may speak of the viscosity of the 
earth but it is not in general the viscosity of a 
fluid. We may speak of internal friction but the 
molecules are bound; they are not free as in a 
fluid and, at least in the upper crust of the earth, 
they are for the most part arranged in crystal 
lattices. Plastic flow probably takes place by 
gliding on cleavage planes and by recrystalliza- 
tion in the presence of more than one phase of 
the chemical system. Thus hydrodynamics when 
applied to the solid earth is somewhat of a 
misnomer. Tectonophysics requires the building 
up of a new mathematical science which we might 
call plastodynamics. The engineers have made a 
good start for us in studying the plasticity of 
metals. The rocks of the earth’s crust as we 
meet them exposed at the surface or as laboratory 
specimens are not only elastically rigid but they 
‘exert a resistance to change of shape beyond 
their elastic limit. They have strength. We may 
define the strength of a rock in different ways, as 
the force required to start rupture, and as the 
magnitude of the largest internal stress system 
which can be set up without plastic flow. We 
shall find that the strength is characteristic of 
the material in a given set of conditions only. 
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We shall find that the term stress cannot be 
applied univocally to the earth and its constit- 
uents. Furthermore, the distinction between the 
terms set, creep, and plastic flow as applied to 
the yielding of a metal, rock, or mineral would 
seem to be arbitrary and the behavior described 
by these terms would seem to differ widely 
from substance to substance. It is not surprising, 
therefore, to find geologists, geophysicists, and 
tectonophysicists divided into rival camps, one 
of which would deny any strength to the rocks 
at considerable depth in the earth, another of 
which would maintain that rocks at all depths 
retain their strength and even increase it greatly 
under pressure, but that the large scale of earth 
phenomena dwarfs the influence of molecule on 
molecule which constitutes strength and inte- 
grates the tiny margin of yielding into so long a 
series that structures composed of the strongest 
rocks behave in the gross like mud, and that the 
earth, though composed of strong materials in 
detail, becomes as a whole so plastic as easily to 
satisfy all requirements of mobility, diastro- 
phism, and orogeny. A third group would hold 
that rocks at the surface have all of the strength 
observed in the laboratory but lose it at the 
depth of a few tens of kilometers largely because 
of high temperature. 

It is nothing new to the meteorologist to say 
that the earth is a thermodynamic machine. 
Neither is it new to the petrologist that thermo- 
dynamic principles are fundamental in the study 
of igneous rocks. Convection currents were a 
matter of course in the old conception of the 
earth as a ball of molten lava with a thin crust. 
But it is quite a different matter to speak of 
thermodynamic systems in an earth that on the 
average is at least as rigid as steel in the labora- 
tory and especially to speak of thermodynamic 
systems in the shell or mantle which lies between 
the crust and the core of the earth and which 
possesses a rigidity or viscosity more than twice 
that of steel. Yet the setting up of large-scale 
convection currents is precisely the explanation 
that is offered by many for the great ocean deeps 
which are so active as earthquake sources. It is 
clear that thermodynamics of fluids cannot be 
applied without modification to such conditions. 
The energy relations in the case of the solid 
earth will be different. 
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In conclusion we may say that the new science 
of tectonophysics is making extraordinarily 
valuable contributions of fact to our knowledge 
of the materials in the crust of the earth. But 
this added knowledge instead of simplifying the 
problem of the origin of earthquakes is immensely 
complicating it. We need facts and more facts; 
but we shall be overwhelmed by their very 


variety unless we also give our attention to the 
development of new and more powerful mathe- 
matical tools which will help us to coordinate 
these facts along the lines of a more general 
theory of elasticity, thermodynamics, and hydro- 
dynamics. It is only by means of such analysis 
that we can see where we are going and which are 
the wider gaps in our knowledge of geodynamics. 


REFERENCES 


(1) L. H. Barnett, “‘The chemistry of the earth’s core,” 
J. Geol. 32, 615-635 (1924). G. Tammann, “Zur 
Analyse des Erdinneren. I. Erdkérper: Zusammen- 
setzung seiner verschiedenen Schichten,”’ Zeits. f. 
anorg. allgem. Chemie 131, 96-100 (1923); ‘‘IT. Auf- 
bau des Erdinneres aus verschiedenen Schichten,”’ 
ibid. 134, 269-296 (1924). V. M. Goldschmidt, 
“Probleme und Methoden der Geochemie,”’ Beitr. 
z. Geophys. 15, 38-50 (1926). , 

(2) Bailey Willis, ““The mechanics of Appalachian struc- 
ture,” Thirteenth Annual Report, U. S. Geological 
Survey, Part 2 (1891-1892) pp. 217-281. M. King 
Hubbert, “‘Theory of scale models as applied to the 
study of geologic structures,’” Bull. Geol. Soc. Am. 
48, 1459-1520 (1937) (contains copious references). 

(3) J. B. Macelwane, “Evidence on the interior of the 
earth derived from seismic sources.’”’ Chapter X. 
Interior of the Earth (McGraw-Hill, 1939). See 
especially Table 25, pp. 223-224; Table 26, p. 225; 
Table 27, p. 226; Tables 28 and 29, p. 227. 

(4) F. D. Adams and E. G. Coker, An Investigation into 
the Elastic Constants of Rocks, More Especially with 
Reference to Cubic Compressibility (Carnegie Institu- 
tion of Washington, Washington, D. C., 1906). 

(5) L. H. Adams and E. D. Williamson, ‘‘The compressi- 
bility of minerals and rocks at high pressures,” 
J. Frank. Inst. 195, 475-529 (1923). 

(6) P. W. Bridgman, The Physics of High Pressure (Mac- 
millan, 1931); ‘‘The compressibility of several 
artificial and natural glasses,” Am. J. Sci. 10, 
359-367 (1925); ‘‘Linear compressibility of fourteen 


VOLUME 14, MARCH, 1943 


natural crystals,’ ibid. 482-498; ‘‘Polymorphism, 
principally of the elements, up to 50,000 kg/cm?,”’ 
Phys. Rev. 48, 893-906 (1935); ‘Shearing phe- 
nomena at high pressure of possible importance for 
geology,” J. Geol. 44, 653-669 (1936); ‘‘The high- 
pressure behavior of miscellaneous minerals,’”’ Am. 
J. Sci. 237, 7-18 (1939). 

Griggs, “Deformation of single calcite crystals 
under high confining pressures,’’ Am. Mineral. 23, 
28-33 (1938) ; ‘Creep of rocks,” J. Geol. 47, 225-251 
(1939); ‘Experimental flow of rocks under condi- 
tions favoring recrystallization,’’ Bull. Geol. Soc. 
Am. 51, 1001-1034 (1940). 

(8) F. Birch and R. R. Law, ‘‘ Measurement of compressi- 
bility at high pressures and high temperatures,” 
Bull. Geol. Soc. Am. 46, 1219-1250 (1935); F. 
Birch and R. B. Dow, ‘‘Compressibility of rocks 
and glasses at high temperatures and pressures: 
Seismological Application,”’ Bull. Geol. Soc. Am. 
47, 1235-1256 (1936); F. Birch, “The effect of 
pressure on the modulus of rigidity of several 
metals and glasses,’ J. App. Phys. 8, 129-133 
(1937); F. Birch and D. Bancroft, “The effect of 
pressure on the rigidity of rocks,” J. Geol. 46, 
59-141 (1938); F. Birch and D. Bancroft, ‘‘Elastic- 
ity and internal friction in a long column of granite,”’ 
Bull. Seismol. Soc. Am. 28, 243-254 (1938); F. 
Birch and D. Bancroft, ‘‘New measurements of the 
rigidity of rocks at high pressures,” J. Geol. 48, 
752-766 (1940). 


(7) D. 


133 














The Interior of the Earth 
Viewed in Relation to Earthquake Causes 


VI. What We Learn from Thermodynamics 


Summary 
By JOSEPH LYNCH, S.J. 


Fordham University, New York, New York 


N closing this series of papers, it may be de- 


sirable to say a word or two about a viewpoint ° 


that has been only indirectly touched on, namely, 
the thermodynamic viewpoint. The only facts on 
which we can base any thermodynamic conclu- 
sions with regard to the nature of the earth’s 
interior are: (1) The existence of a thermal gra- 
dient amounting to about 30 degrees per kilom- 
eter. We have proof only of the existence of such 
a gradient to the depth of some few kilometers. 
The presumption is that the gradient continues 
below this depth but it is only presumption. It 
seems, however, a reasonable assumption that 
the gradient should exist all the way to the center 
but with decreasing steepness. This thermal gra- 
dient multiplied by the conductivity gives us the 
outward surface heat flux. (2) The Clapeyron- 
Clausius equation tells us that the melting point 
of rocks, averaging about 1400 degrees C at the 
earth’s surface, increases with pressure and hence 
with depth. The increase amounts to about 3 
degrees C per kilometer. Here again we have no 
assurance that such increase continues through- 
out the earth, but Bridgman on the basis of his 
experimental work with high pressures states that 
the indications are that it does so continue but 
probably at a slower rate. These are the only 
facts we have to go on. If, however, we make the 
assumption that the earth was originally a mol- 
ten mass, then as shown by Jeffreys, the tem- 
perature gradient of such a liquid sphere, cooling 
adiabatically by convection, would be about 
three-tenths of a degree per kilometer. The melt- 
ing point gradient is ten times greater than this, 
namely, 3 degrees rise per kilometer. Hence, 
solidification would start not at the surface but 
at some depth below the surface. The precise 
depth would be determined by the initial tem- 
perature of the molten sphere. Jeffreys and others 
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now believe that solidification began halfway 
down in the earth at what is generally accepted 
as the core of the earth. 

Because of their high density, most radioactive 
materials must have been imprisoned in this 
liquid core when solidification began and as Hess 
pointed out in his paper, it is hard to understand 
how they could all have escaped, as Jeffreys and 
Holmes suggest, in the form of volatile matter; 
hence, one would expect some radioactive ma- 
terial still to be in the core. 

If we take the age of the earth as about two 
billion years, Slichter has shown in the paper 
already mentioned by Hess, that thermal phe- 
nomena in the earth are practically at a stand- 
still. He has shown, moreover, that a layered 
earth compatible with the findings of seismology 
and with radioactivity, continuing with decreas- 
ing intensity down to the center, would have a 
remarkably low temperature. Depending on the 
values chosen for conductivity and diffusivity, 
reasonable figures for the temperature of the core 
would range from somewhat below to somewhat 
above 2000 degrees C. Assuming a rise in the 
melting point of rocks of 3 degrees C per km, the 
solidification temperature, assuming 1400 degrees 
for the surface, would be 2000 degrees C at a 
depth of only 200 km, so that at the core it should 
be higher than 2000 degrees C and the core con- 
sequently solid. Tammann would probably object 
to this that the melting point, after increasing at 
first with depth, actually decreases with further 
depth instead of increasing. To this objection 
we can only repeat the statement of Bridgman 
that there is no experimental evidence to support 
Tammann’s view—on the contrary what evi- 
dence there is contradicts it. Thermodynamics 
would seem, therefore, to argue for an earth solid 
right through to the center. 
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Such a picture of the interior of the earth 
would fit in better than any other with the earth’s 
magnetic field. As Lynch has shown, it is im- 
possible to ascribe the magnetic field of the earth 
to any mere ferromagnetic shell. As mentioned 
in the same paper, a core temperature of 2200 
degrees C according to McNish would probably 
permit of a ferromagnetic core predominately 
nickel. Such a solid ferromagnetic core, would be 
more acceptable than a liquid core to the ter- 
restrial magnetician. 

Hess has shown a big discrepancy between the 
low radioactive figures of Evans and Goodman 
and the known value of gamma-radiation over 
rocks. He inclines to the higher values of Holmes 
and Jeffreys pending further investigation, and 
prefers as the best picture of the earth’s interior 
so far suggested a layered earth, with radio- 
activity decreasing discontinuously in the layers 
all the way down through the core, the tempera- 
ture of the latter being around 2000 degrees C. 
Thus, thermodynamics, terrestrial magnetism, 
and radioactivity seem to favor an earth solid 
right through to the center. 

Seismology certainly demands a solid mantle 
since both compression and shear waves traverse 
it with ease. Hodgson in his paper on seismology 
rightly confined himself to a discussion of the 
mantle and omitted the core. As he pointed out, 
no earthquakes are as yet known to have oc- 
curred at a greater depth than 700 km and hence 
the part of the earth’s interior most connected 
with earthquake cause would be this upper part 
of the mantle. However, a very natural question 
would be—does anything in seismology preclude 
the possibility of a quake as deep seated as the 
core itself? Ultimately this resolves itself into 
whether or not seismology precludes the possi- 
bility of a liquid core. The only seismic tool to 
throw light on that question is the shear wave. 
Whether or not the shear wave traverses the core 
is still a debated question. Bastings in New 
Zealand after exhaustive study is convinced that 
it does. Krumbach in Germany, Imamura in 
Japan, and Macelwane in America incline to the 
same view but with slightly less conviction. All 
are agreed that the shear wave is robbed of most 
of its energy in so passing or attempting to pass 
through the core, but the evidence from seis- 
mology seems to favor a solid rather than a liquid 
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core. Geodesy and tectonophysics throw no light 
directly on the core, their findings being equally 
well satisfied by a solid or a liquid core. Geodesy, 
however, in as much as it postulates isostasy, 
requires a region of flow in which compensation 
can be effected. According to Jeffreys such plastic 
flow is chiefly the result of relaxation resulting 
from temperature molecular vibration. Since this 
is a function of temperature and necessarily in- 
creases with temperature, one would expect on 
the basis of the known temperature gradient of 
30 degrees per km that conditions for flow would 
increase as we descend into the earth. Why ideal 
conditions for flow should exist just at the iso- 
static layer is not indicated in any of the papers 
in the present series, but Heck’s paper shows that 
geodetic findings are best satisfied by assuming 
a compensating layer about 60 miles down. To 
the man in the street it seems paradoxical to say 
that we can have plastic flow in a mantle which 
is yet rigid enough to transmit shear waves. 
Jeffreys’ answer to the paradox is that time inter- 
vals such as occur in the passage of a shear wave 
are too short to permit of relaxation coming into 
play—hence, for rapid phenomena the material 
shows rigidity. The much longer time interval 
involved in rock flow, however, permits relaxa- 
tion and consequently flow. 

From the topics discussed in this series of 
papers, it is apparent that the nature of the 
earth’s interior is still in doubt. The evidence 
seems to point toward a solid core. Yet tidal 
phenomena and the Eulerian period call for a core 
of low rigidity, so that we need a solid earth with 
a solid core of low rigidity. Two years ago the 
writer advanced as the possible solution of such 
a rigidless solid core, a solid solution—experi- 
ments with palladium showing beyond doubt 
that the effect of absorbed hydrogen in the pal- 
ladium was to lower the rigidity of the metal. 

With regard to the cause of earthquakes, as 
Heck pointed out, we must distinguish between 
the proximate and remote causes. Both Heck and 
Hodgson agreed that the immediate or proximate 
cause of earthquakes is slippage along a fault 
plane. Several cogent reasons were advanced by 
Hodgson to indicate that deep focus quakes do 
not materially differ from shallow quakes in this 
immediate mechanism. In both papers it was 
agreed that the best cause so far suggested for 
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this slippage was the elastic rebound theory of 
Reid, fair confirmation of which is had from 
geodesy. As to the more remote causes that set 
up the stresses which produce this rebound, it 
was agreed that beyond the fact that they are 


be said of them. Contraction, erosion, and deposi- 
tion, with resulting isostatic compensation were 
agreed on as probably playing a prominent part 
in earthquake production, but the complete cause 
is a jigsaw puzzle which we have not yet been 





connected with mountain building, little could able to put together. 


136 





NE of the influences of the war upon physicists is to make them conscious 

of applied physics as an integral part of their interest. Using the term phys- 
ics in this more comprehensive sense, one can easily recognize that there must 
be a connection between the recent advances of physics and the structure of our 
social and political institutions of tomorrow. Physics will prolong life. It will 
give greater physical comfort to human beings. It will make available more 
power per capita. It will enable us to hear more, see more, travel farther, and 
increase the possible experiences of the individual in his lifetime. It will increase 
the contact between nations, thus making tradition and past history of less 
and less dominating importance in nationalism. Even though the spirit of in- 
tense nationalism has increased during the period of industrial development, 
advances in physics and its applications will have a definite tendency to reduce 
extreme nationalism and to facilitate the adjustment of human relations. 
Physics, along with the other sciences, gives confidence in the use of the human 
mind and in the mind’s creativeness. We need this confidence to offset the 
backward view of those who imagine we are rushing toward destruction.—G. 
W. STEwart, in “Social Implications of Physics,” in the American Journal of 
Physics, February, 1943. 
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Résumé of Recent Research 








High Temperature 
X-Ray Camera* 


Some organic com- 
pounds exhibit phase 
changes with respect to 
temperature as shown by the x-ray diffraction 
camera.' Sodium stearate is such a compound. 

The sample is mounted in a thin glass capillary 
not over 1 mm in diameter and inserted in the 
sample holder designated as (1) in Fig. 1. The 
sample holder shaft is of stainless steel to mini- 
mize conduction. The sample holder is made of 
brass. A heater coil (2) wound on sheet asbestos 
on a stainless steel cylinder is above and below 
the sample to help decrease sharp temperature 
gradients. A small glass wool insulated thermo- 
couple placed not over 3 mm from the sample 
serves for temperature measurement. Easy access 
can be had to the oven coils. A water-cooled 
jacket (3) surrounding the oven prevents hot air 
transfer to the x-ray film on the circumference 
of the camera. The camera takes about 20 min- 
utes to reach thermal equilibrium. 

The filtered Cu K x-ray beam enters the water 
jacket from the pinhole system (4) and is 
diffracted by the sample. The exit window (5) in 
the water jacket permits passage of the forward 
diffracted beam 90° each side of the direct x-ray 
beam. The exit window is covered with thin 

















* Alexander de Bretteville, Jr., “‘An x-ray high tempera- 
ture camera,”’ Rev. Sci. Inst. 13, 481 (1942). 

1A. de Bretteville, Jr.. and J. W. McBain, ‘X-ray 
diffraction investigation of sodium stearate from room 
temperature to the melting point,’’ to be published. 
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Cellophane fastened with Scotch tape except for 
an 8-mm hole for the direct beam. The Cellophane 
prevents the escape of heat and thus furnishes 
additional protection for the x-ray film. 

The camera body is essentially a cylinder of 
57.3-mm radius into which the cover unit, con- 
taining the water jacket, heater coil, and sample 
fits concentrically. The cover unit is fastened 
down around the circumference of the camera 
body by screws. 

The film is centered on the circumference of 
the camera with respect to the direct x-ray beam 
by a tapered aluminum plug inserted in the exit 
port. This plug is then removed and the exit port 
cap is replaced. 

It is estimated that a temperature as high as 
500°C could be obtained. 





New Instrument Booklets 








Bakelite Review, July 1942, Vol. 14, No. 2, describes the 
presentation of the Navy ‘‘E”’ awards to the Bound Brook 
plant of Bakelite Corporation and Wyandotte, Michigan 
plant of Halowax Corporation. Interesting articles on the 
following topics are featured: ‘Plastic parts for machine 
tools,’ ‘Finishing of resin-bonded plywood for aircraft 
construction,”’ “Heating with high frequency,” ‘Special 
molds and molding technique essential for fabrication of 
high impact-resistant plastics.’’ Quarterly publication of 
Bakelite Corporation, Unit of Union Carbide and Carbon 
Corporation, 30 East 42nd Street, New York, New York. 


Cenco News Chats, September, 1942, No. 38, in its usual 
informal style presents information about the construction 
and performance of selected Cenco products. The section 
entitled ‘‘Cencograph”’ contains views of the buildings and 
research activities of Parke, Davis & Company. An edi- 
torial, ‘Appraising science anew," tresses the peacetime 
contributions and responsibility of science. Occasional 
publication of Central Scientific Company, 1700 Irving 
Park Road, Chicago, Illinois. 


Chemical Specialties presents a list of reagents most in 
demand for analytical work in industrial, hospital, and 
educational laboratories. Each reagent bears a printed label 
which states the exact percentage of impurities in the lot; 
the reagents meet with American Chemical Society 
specifications. This booklet and the products listed in it 
are available from: Eimer and Amend, 635 Greenwich 
Street, New York, New York; Fisher Scientific Company, 
717 Forbes Street, Pittsburgh, Pennsylvania; Fisher 
Scientific Company, Ltd., 904 St. James Street, Montreal, 
Quebec, Canada (95 pages). 
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Effective Use and Proper Care of the Microtome written 
by Oscar W. Richards, 1942, is thoroughly illustrated. The 
author summarizes the available facts on microtomes and 


describes the use and care of the microtome. He describes 
special experimental work which he did on section cutting. 
Copies may be obtained at a cost of 25 cents by writing to 
Spencer Lens Company, Buffalo, New York (88 pages). 


General Radio Experimenter, Vol. 16, No. 6, November, 
1942, discusses measuring instruments which allow accu- 
rate evaluation of the various factors involved in receiver 
performance. A new reasonably priced generator called the 
Type 805-A Standard-Signal Generator is available; it is 
designed to satisfy the need for an amplitude modulated 
signal generator in the frequency range up to 50 megacycles 
that would be suitable for testing all types of receivers 
ranging from special military equipment to high-fidelity 
broadcast sets. Published by General Radio Company, 30 
State Street, Cambridge A, Massachusetts (5 pages). 


Harder Tin-Base Alloys is a new technical publication 
on methods of obtaining new and stronger tin-base alloys 
suitable for use as bearing metals. The booklet may be 
obtained from: Tin Research Institute, Fraser Road, 
Greenford, Middlesex, England. The Tin Research Insti- 
tute is controlled by the International Tin Research and 
Development Council. 


Inco, Vol. 18, No. 3, Summer-Fall edition, presents 
articles on various services offered by the International 
Nickel Company and describes new and outstanding uses of 
Inconel, monel, and other nickel alloys. Publication of the 
International Nickel Company, Inc., 67 Wall Street, New 
York, New York (32 pages). 


Process Industries Quarterly, Vol. 7, No. 2, 1942, 
emphasizes ways to overcome wartime shortages. Since 
existing equipment must be made to last for the duration, 
it is suggested that sheet nickel linings be used on steel 
units. There is an interesting article on reconditioning worn 
pump shafts and rods. Publication of the International 
Nickel Company, Inc., 67 Wall Street, New York, New 
York. 





New Books 


A.S.T.M. New Standard Specifications 


Through action of its Committee E-10 on Standards, the 
American Society for Testing Materials has recently issued 
some 32 new standard specifications and test methods. Of 
these, 13 are emergency specifications and tests designed to 
expedite procurement or conservation of critical or strategic 
materials. Emergency Specification ES-20 covers malleable 








iron flanges, pipe fittings, and valve parts including parts to . 


be assembled, manufactured in advance, and supplied from 
stock by the manufacturer or distributor. Twelve of the 
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emergency specifications cover various types of thermal 
insulating material in different forms and for service at 
widely varying temperatures ranging from surface tempera- 
ture of 200 to 1900°F depending on type of materials. The 
serial designations and brief titles are as follows: ES-8, 85 
percent magnesia cement; ES-9, long fiber asbestos cement; 
ES-10, mineral wool thermal cement; ES-11, expanded or 
exfoliated mica cement; ES-12, diatomaceous earth cement 
(600 to 1200°F); ES-13, diatomaceous earth cement (1200 
to 1900°F); ES-14, blanket insulating for building purposes; 
ES-15, blanket insulating for industrial purposes; ES-16, 
blanket insulating for refrigeration; ES-17, preformed pipe 
covering; ES-18, preformed block insulating; ES-19, 
structural board insulation. 

The Society also announced that approval had been 
given to modification in the emergency specifications for 
carbon-chromium ball and roller bearing steels, ES-5. 

Eighteen of the new tentative specifications and tests 
announced were developed in regular course by the various 
A.S.T.M. technical committees and deal with: recom- 
mended practice for making torsion tests of cast iron; 
permeability of feebly magnetic materials, permeability 
and core loss of flat-rolled magnetic materials at low 
alternating inductions using 28-cm specimen, and perme- 
ability and core loss of flat-rolled magnetic materials using 
28-cm specimen; new tentative methods for evaluating 
paints; tests for compressive strength, coefficient of linear 
thermal expansion, and water vapor permeability of 
plastics; specifications covering gravel for bituminous con- 
crete base and surface courses and for water-bound base and 
surface courses. 

All of these new specifications are being printed in 
separate pamphlet form (the emergency specifications on 
pink stock) and copies can be obtained at 25 cents each 
from A.S.T.M. Headquarters, 260 S. Broad Street, Phila- 
delphia, Pennsylvania. They will also be bound in the 1942 
A.S.T.M. Book of Standards. 





Here and There 


Enemy-Owned Patents Available for Industry 








The United States Government, through the Office of 
Alien Property Custodian, now holds about 50,000 patents 
formerly owned by residents of enemy and enemy-occupied 
countries. The inventions disclosed in these patents repre- 
sent research achievements of great potential value. The 
patents are being offered for use by American industry to 
aid in winning the war and as a contribution toward 
rebuilding our peace-time economy. 

A pamphlet has been prepared to provide a brief outline 
of the policies which have been adopted for the adminis- 
tration of these patent holdings. There is also included an 
index to the classified lists of vested patents for use as a 
guide in ordering specific sections. A copy may be obtained 
by writing directly to the Office of Alien Property Cus- 
todian, Washington, D. C. 
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Honors and Awards 


William L. Batt, Chairman of the Engineering Industrial 
Division of the National Research Council, was awarded an 
honorary degree of Doctor of Sciences by the Drexel Insti- 
tute of Technology at Philadelphia on December 17, 1942. 


Sir Robert Robinson was awarded the Copley Medal of 
the Royal Society for his outstanding work in organic 
chemistry. 
Army-Navy “E” Award was presented by Major 
General Benjamin W..Chidlaw, Assistant Chief of Staff of 
the Engineering and Materials Division of the United 
States Army Air Forces, to the Nylon Research Laboratory 
and Pilot Plant of the du Pont Company on December 29th. 


* 


Bequest to Northwestern University 


A bequest of more than $20,000,000 was made to 
Northwestern University by the will of Walter P. Murphy, 
manufacturer of railway supplies who died on November 
16, 1942. The fund is to be used to develop, maintain, and 
operate the 
University. 


Technological Institute. of Northwestern 


* 


Joint Council of Professional Scientists 


A Joint Council of Professional Scientists has been 
formed in England to aid in the training and use of scien- 
tists in the war effort. Some of the members of the council 
are: representing physicists; Sir Lawrence Bragg, Professor 
J. A. Crowther, Mr. E. R. Davies, Dr. B. A. Keen, Dr. H. 
R. Lang, Dr. H. Lowery; representing mathematicians, 
Professor S. Chapman; representing geologists, Professor 
H. H. Read. 


* 


File on Foreign Scientific and Technical Publications 


Harold Lancour, Librarian of Cooper Union, Chairman 
of Engineering School Libraries Section of the Association 
of College and Reference Libraries, reports that the Library 
of Congress is now attempting to make a master file of 
current technical and scientific literature published abroad 
and obtainable in this country. This file is expected to aid 
many phases of war research. 


* 


Staff Changes 


D. L. Beeman has been appointed engineer in the 
industrial power section of General Electric Company 
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succeeding E. G. Merrick and R. S. Sage has been made 
engineer in the mining section succeeding F. L. Stone. 


Fitzhugh W. Boggs is now research engineer at the 
laboratories of Westinghouse Electric and Manufacturing 
Company, East Pittsburgh. Elizabeth M. Boggs is a 
lecturer in the department of physics, at the University of 
Pittsburgh. 


Harold W. Coles of the Mellon Institute has become 
associated with the research staff of the Bausch & Lomb 
Optical Company, Rochester, New York. He has been 
placed in charge of the new department of organic 
chemistry. 


Wilmer T. Rinehart, recently associated with the United 
States Gypsum Company, has joined the chemical engi- 
neering staff of the Armour Research Foundation. 


Sir John Russell, known for his work in agricultural 
science, will retire October 30th from the directorship of the 
Rothamsted Experimental Station. The Station reaches its 
100th year in 1943. 


William W. Skinner has been appointed chief of the 
Bureau of Agricultural Chemistry and Engineering of the 
United States Department of Agriculture, to fill the 
vacancy left by the death of Henry G. Knight. 


* 


Technical Personnel Needed 


The Naval Ordnance Laboratory at the Navy Yard, 
Washington, D. C., is a research and development agency 
of the Bureau of Ordnance. This laboratory needs addi- 
tional technical personnel. Physicists and electrical engi- 
neers with electronics experience are required; needed also 
are mechanical engineers familiar with the design of small 
mechanical movements or mechanisms. Several openings 
are available in technical report writing and editing. 


* 


Bulletin of Mathematical Biophysics 


The March, 1943 issue of The Bulletin of Mathematical 
Biophysics, edited by N. Rashevsky and published by the 
University of Chicago Press, will contain the following 
articles: 


Some Theoretical Consideration Concerning the Interchange of 
Metabolites between Capillaries and Tissue. INGRAM BLocn. 

Theory of Cellular Pulsations. HENry E. STANTON. 

The Linear Theory of Neuron Networks: The Dynamic Problem. 
WALTER PIrTTs. 

Outline of a New Mathematical Approach to General Biology: | 
N. RASHEVSKY. 
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Necrology 


Frederick M. Beckett, consultant to the Union Carbide 
and Carbon Corporation died December 1 at the age of 67. 


Purnendu Nath Chakravorty, native of India, research 
chemist with the Upjohn Company and formerly associated 
with the Department of Chemistry at Princeton Uni- 
versity, died December 23 at the age of 37. 


Clyde Coleman, research chemist for the United States 
Rubber Company, died November 6 at the age of 59. 


Alfred Nelson Finn, for many years chief of the depart- 
ment of optical glass of the National Bureau of Standards, 
died at the age of 60. 


Dr. Arthur P. Honess, Professor of Mineralogy and 


Petrology at The Pennsylvania State College, died October 
17 at the age of fifty-five. 
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James Edmund Ives, senior physicist of the United 
States Public Health Service until his retirement in 1936, 
died January 1 at the age of 70. 


Andrew H. Palmer, meteorologist for the United States 
Weather Bureau from 1914 to 1924, and later Superin- 
tendent of Crops at the Aetna Affiliated Companies, San 
Francisco, died December 26 at the age of 56. 


Captain George Tyrrell McCaw, authority on Geodosy, 
Editor of the Empire Survey Review, and _ technical 
assistant to the Geographical Section of the General Staff, 
War Office of England, died October 17 at the age of 72. 


Robert Peele, Professor emeritus of Mining Engineering 
of the School of Mines of Columbia University and Editor 
of The Mining Engineer's Handbook since 1917, died De- 
cember 8 at the age of eighty-four. 


Herman Stabler, chief of the conservation branch of the 
United States Geological Survey since 1925, died on No- 
vember 24 at the age of sixty-three. 
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The Lubrication by Thin Metallic Films and the Action of Bearing Metals 


F. P. BowpDEN AND D. TABOR 
Council for Scientific and Industrial Research, East Melbourne, Australia 


(Received November 21, 1942) 


An investigation has been made of the role of thin metallic films in reducing the friction 
and wear between metal surfaces. It is shown that with metallic films possessing suitable me- 
chanical and surface properties the coefficient of friction is very low and may be similar to 
that observed on ice. In many cases the behavior of these thin metallic films resembles that 
of ordinary lubricants except that Amontons’ law is not obeyed. The friction and wear prop- 
erties of some copper-lead alloys have been studied and the results indicate that the lead in 
the alloy is extruded during sliding and forms a thin lubricating film on the hard copper 
matrix. It is suggested that these metallic films play an important part in the action of many 


bearing alloys. 


INTRODUCTION 


HE last paper' discussed the role of shearing 

and ploughing in the friction of sliding 
metals and described the lubricating properties 
of thin indium films. This work has been ex- 
tended to films of lead and copper. The part these 
films may play in the action of certain bearing 
metals is also considered. It has been maintained 
in the past that an essential characteristic of a 
bearing alloy is that it should possess a duplex 
structure consisting of hard crystals embedded 
in a relatively soft matrix. It is suggested (see, 
for example, Desch*) that the function of the 
hard crystals is to resist wear and that of the 
softer constituents to permit a more uniform 
distribution of the load, by allowing any of the 
hard crystals which are heavily loaded to sink 
so that the load is spread over a greater area. It 
is also suggested that the hollows worn in the 
softer material serve as reservoirs for the lubri- 
cating oil. There is probably something in these 
suggestions and it is certainly true that many 
successful bearing alloys do possess a structure 
of this type. It is clear, however, that it is not a 





1 Bowden, Moore, and Tabor, J. App. Phys. 14, 80 
(1943). 


2 Desch, Metallography (Longmans, Green and Company, 
New York, 1937). 
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complete explanation and in many cases is cer- 
tainly not true at all. With many modern bearing 
alloys, the surface layer does not possess a duplex 
structure at all, but consists of a single pure 
metal. Another very wide class of bearing alloys 
consists of a continuous matrix of the harder 
metal with a small amount of the softer metal 
finely dispersed through it. Some of the copper 
lead alloys provide a good instance of this (see, 
for example, the copper lead alloys shown in Figs. 
10, 11). In this case, the hard copper forms the 
continuous phase so that it is not possible for it 
to “sink’’ into the lead. It is suggested that in 
bearing alloys of this type a thin film of the 
softer metal is extruded by the sliding process 
and acts as a lubricant film on the surface of the 
harder matrix. 


LUBRICATION BY THIN METALLIC FILMS 
The Relation Between Track Width and Friction 


In these experiments, films of indium, lead, 
and copper of varying thicknesses were plated on 
to flat metallic surfaces. The substrate was steel, 
nickel, copper, or silver. The upper surface was 
a hemispherical steel slider of radius of curvature 
0.3 cm. It was found that the friction is almost 
independent of the underlying metal and of the 
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Fic. 1. Friction of steel slider on lead films deposited on 
steel, silver, and copper substrates, as a function of track 
width. It is seen that the friction is determined primarily 
by the track width, whether the friction is determined on 
the surface film or on the bulk lead. 


load. The main effect of increasing the load is to 
cause a slight increase in the deformation of the 
underlying metal, and so cause a slight increase 
in the track width d. Since silver is softer than 
steel the same load will produce a greater defor- 
mation. The track width d and hence the fric- 
tional force F for any given load and film thick- 
ness will, therefore, be greater on silver than on 
steel. This effect was observed. In Fig. 1, F is 
plotted against d for lead films deposited on 
various substrates; the different values of the 
track width being obtained by varying the 
thickness of the deposited layer and the radius 
of curvature of the upper surface. It will be seen 
that in all cases the frictional resistance is deter- 
mined, primarily, by the track width, whether 
the friction is measured on the surface film or on 
the bulk metal. Essentially similar results were 
obtained for indium and copper films. 


The Shear Strength of the Metal Junction 


The real area of contact A is proportional to 
the square of the track width d. It is found, as 
has already been pointed out in the first paper, 
that when F is plotted against d? the results lie 
approximately on a straight line, whether the 
track is made by sliding on the surface film or on 
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the bulk metal. We should, therefore, expect 
that the friction on the different metallic films 
would be proportional to the shear strength of 
the metal. It will be seen from Fig. 2 for indium, 
lead, and copper films that this is so. For any 
given track width the friction of copper>lead 
> indium. 

The values of the shear strengths are given in 
Table I. The second column shows the values 


TABLE I. 


Shear strength g/mm+? ; 
_. Calculated from From shearing of 
friction measurements pure metals 











Indium 325 220 
Lead 1,600 750 
Copper 28,000 16,000 
Steel ball 140,000 90,000 


calculated from the friction measurements (neg- 
lecting the ploughing term P). The third column 
shows the shear strength obtained by the direct 
experimental method of shearing small cylinders 
of the pure metal. For completeness, the values 
obtained for a steel ball sliding on hard steel are 
included at the bottom of the table. 

The proportionality is at once apparent. It 
will be noted that the shear strength of the junc- 
tions calculated from the frictional measurements 
is of the same order of magnitude as that of the 
pure metals but is somewhat greater. This is 
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__ Fic. 2. Friction of steel slider on indium, lead, and copper 
films, as a function of track width. It is seen that for any 
given track width the friction of copper >lead >indium. 
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Fic. 3. Frictional behavior of indium film on steel as wear proceeds. 


probably due to the fact that there was an 
appreciable accumulation of metal in front of the 
slider so that there is some uncertainty in the 
value of A, and to the work-hardening and 
pressure effects described earlier. 


Breakdown and Wear of Films 


If very heavy loads are used, the values of F 
deviate from the curves shown in Fig. 1. After a 
certain load has been exceeded the friction 
increases markedly with load, even though the 
track width remains essentially unaltered. An 
examination of the track shows that this is con- 
nected with the progressive breakdown of the 
film, as the heavier loads are applied. The load 
at which this breakdown occurs depends upon 
the film thickness and upon its strength of ad- 
hesion to the underlying surface. It is also 
influenced by the shape of the slider and takes 
place more readily if this has a small radius of 
_ curvature. 

The film may also be worn away if the slider 
traverses the same track a sufficient number of 
times. Some results obtained with an indium 
film 4X10-* cm thick deposited on tool steel 
(hardness 800 V.D.H.) are shown in Figs. 3 and 4. 
The frictional behavior is shown in Fig. 3 and it 
is seen that initially the motion was smooth and 
the coefficient of friction was about »=0.08. 
With successive slidings over the same track the 
friction gradually rose and after the seventh run 
small stick-slips set in. The friction and the size 
of the stick-slips increased and, after 20 runs the 
friction had risen to n.=0.4. The corresponding 
tracks are shown in Figs. 4 a, b, c, d. It will be 
seen that after the seventh run the steel surface 
is partially exposed and the original lapping 
marks can be seen. After the twentieth run there 
is considerable breakdown of the film. 
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Similar results were obtained with indium 
films deposited on silver surfaces. With a film 
10 cm thick the friction rose with successive 
slidings, and stick-slips set in after the fourth 
run. With a film three times as thick, the friction 
rose more gradually and stick-slips did not begin 
till the eighteenth run. 

The rise in friction with wear of indium films 
is shown in Fig. 5, and it is seen that the thicker 
film is more resistant to wear. This behavior may 
be compared with the dotted curve which repre- 
sents the rate at which stearic acid layers are 
worn off a steel surface. Bowden and Leben have 
shown’ that with 9 molecular layers of the acid 
an appreciable rise of the friction occurs after 
20 successive runs. With 53 molecular layers, 
i.e., about 10-5 cm, no appreciable change was 
observed after 100 successive runs. It is clear 
that the metallic films resist wear quite well, but 
they are more easily worn off the surface than a 
much thinner film of the fatty acid. One factor 
that may be of importance here is the surface 
mobility of the film. The fatty acid molecules 
have a certain mobility and are able to move over 
the surface and repair the damaged film. Many 
metals also possess this ability to wander over 
solid surfaces‘ > but it is probable that the rate 
at which this process occurs is less rapid for solid 
metallic films than it is for the fatty acids. 


The Limiting Film Thickness 


A few experiments were carried out to deter- 
mine the minimum thickness of the deposited 
film that would influence the friction. The results 
for indium on tool steel are shown in Fig. 6. The 


3 Bowden and Leben, Phil. Trans. A239, 1 (1940). 

4 Cockcroft, Proc. Roy. Soc. A119, 293 (1928). 

5 Adam, Physics and Chemistry of Surfaces (Oxford 
University Press, London). 
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(c) 
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Fic. 4. Photomicrographs of track X 34 showing progressive wear of indium film on steel: (a) After first run; (b) After 
second run; (c) After seventh run; (d) After twentieth run. 


film thickness was calculated from the quantity 
of electricity used to deposit it. The upper curved 
slider was of hard steel (radius 0.3 cm) and the 
load was 4000 g. As one might expect, the friction 
decreases as thinner films are used because the 
area of contact A becomes smaller. There is, 
however, a limit to this, and a minimum friction 
is reached when the thickness is of the order of 
10-° cm. With thicknesses less than this, e.g., 
10-* cm (or 50 atomic layers) the film ceases to 
be effective. 

The underlying surface of the tool steel was 
not highly polished nor was it homogeneous. 
These factors might easily cause the indium to 
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be deposited preferentially in certain regions, 
leaving patches of the surface uncovered or only 
very thinly covered. It is possible that if highly 
polished surfaces of a uniform metal were used 
as a substrate, very much thinner layers, even 
perhaps one or two molecular layers, might prove 
effective. 


The Effect of Temperature 


The effect of temperature on the friction of an 
indium film is shown in Fig. 7. The friction 
decreases steadily as the temperature rises, and 
reaches a minimum at the melting point of the 
indium. When the melting is complete there is a 
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sharp rise in friction and, if the film is allowed to 
solidify, the friction again falls to a low value 
(see Fig. 8). 

Similar results were observed with a lead film 
3X 10-* cm thick deposited on to tool steel. The 
friction fell to about half its value over the tem- 
perature range 20° to 270°C. At the melting 
point there was a further slight fall followed by a 
sudden rise as the film became completely liquid. 

The decrease of friction with increasing tem- 
perature is of interest, and should be contrasted 
with the effect of temperature on the friction of 
pure metals. Earlier experiments® showed that 
when the friction was measured in the usual way 
by applying a load to the metals, the effect of 
temperature was relatively small. As was pointed 
out at the time, this is to be expected. As the 
temperature is raised the flow pressure of the 
metal decreases, so that the real area of contact 
A of the metallic junctions will increase. Since 
the shear strength s of the metals falls off propor- 
tionally, however, the force required to break the 
junctions may be approximately the same, and 
the friction may be nearly independent of the 
temperature. In the experiments with the thin 
films, however, the real area of contact A remains 
almost constant when the temperature is raised, 
since the load is taken primarily by the hard 
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Fic. 5. Rise in friction with wear of indium films. For 
comparison, the rate of wear of nine molecular layers of 
stearic acid on steel is shown in the dotted curve. 


6 Bowden and Hughes, Proc. Roy. Soc. A172, 263 (1939). 
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Fic. 6. Effect of film thickness on the lubricating properties 
of thin indium films deposited on tool steel. 


underlying metal. We are therefore, able to 
evaluate separately the effects of A and s and to 
observe the decrease in F as s decreases. 

With films of both lead and indium, F falls to 
about half its value as the temperature is raised 
from room temperature to the melting point. 
The decrease in F is almost linear with tem- 
perature. If we consider the case of the lead film 
and call the frictional force at 20° Feo and at 
300° F300, we may define y the temperature coef- 
ficient of F by the equation 


F300= Foo(1 +t) where t= 300— 20, 
then 
y= —0.0015. 


With the indium films over the temperature 
range 20° to 150°C, y= —0.0025. Very close to 
the melting point when incipient fusion begins 
the friction falls still further, and it is in this 
region that F reaches its lowest value. Once 
melting is complete there is a sharp rise and F 
increases to a high value. In the case of indium 
deposited on tool steel the value reached is almost 
as high as that of uncoated tool steel and it is 
clear that the molten film affords little protection 
to the surface. Since the molten indium does not 
readily wet the tool steel surface, this is not sur- 
prising. With the lead film on the silver surface 
a similar abrupt rise in friction was observed but 
the extent of the rise was very much less. The 
friction of the steel slider on the silver surface in 
the presence of molten lead was only about one 
quarter of that characteristic of steel on clean 
silver. Molten lead readily wets the silver, and 
it is clear that it is still able to afford some pro- 
tection to the surface and to act as a lubricant. 
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Fic. 7. Effect of temperature on the frictional behavior of indium film 4X 10~* cm thick deposited on silver. 
Heating up. 


Some subsidiary experiments carried out on 
mercury films are of some interest in this con- 
nection. [t was found that minute quantities of 
mercury on silver produced a very marked re- 
duction in friction. The mercury readily wets the 
silver and forms an amalgam with it, and this 
apparently enables it to act as a fairly effective 
lubricant. These observations show the important 
part which surface tension effects may play in the 
action of bearing metals, and in the friction and 
lubrication of other combinations of metals. 


Metallic Films as Lubricants 


Some experiments were carried out to deter- 
mine the extent to which the coefficient of friction 
u is independent of the load (Amontons’ law). 
The results for tool steel surfaces are shown in 
Fig. 9 for loads varying from 400 to 8000 g. 
Curve I is for the unlubricated steel; Curve II 
for steel lubricated with a film of mineral oil; 
and Curve III for steel lubricated with a thin 
film of indium approximately 4X10-* cm thick. 
The comparative behavior is striking. The coef- 
ficient of friction for the unlubricated steel, and 
for the steel lubricated with the mineral oil is 
independent of the load. With the indium film, 
however, uw decreases markedly as the load 
increases. 

In the case of unlubricated steel, or of steel 
lubricated with the mineral oil, the area of 
contact should be proportional to the applied 


load, since it depends on the amount of plastic 
deformation that occurs. For this reason F should 
be directly proportional to the load, i.e., Amon- 
tons’ law should hold and the coefficient of 
friction should be constant. This is in fact ob- 
served. For the unlubricated steel surfaces, the 
coefficient of friction is constant and p=0.34, 
while for the surfaces lubricated with the mineral 
oil, the coefficient of friction is again constant 
and w=ca. 0.14. 

If a thin film of a soft metal is used as a lubri- 
cant, however, A will no longer be proportional 
to the load. The real area of metallic contact now 
increases only to a slight extent when the load is 
increased, since the increase in A due to the 
increased deformation of the underlying steel is 
small compared with its actual value A. The 
shear strength s remains constant and is approxi- 
mately that of metallic indium, so that F should 
be almost independent of the load. Amontons’ 
law will no longer hold, and the coefficient of 
friction » will not be constant but will decrease 
as the load increases. It will be seen from Fig. 9 
that this occurs. The coefficient of friction falls 
from n»=0.2 at the light load to ».=0.04 when the 
heaviest load is used. If it were possible to 
lubricate the surfaces with metallic films of 
molecular dimensions, we should expect a closer 
agreement with Amontons’ law. 

In some ways the behavior of the metallic films 
closely resembles that of ordinary lubricant films. 
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Fic. 8. Effect of temperature 
on the frictional behavior of 
indium film 4X10-* cm_ thick 
deposited on silver. Cooling 
down. 














146 


JOURNAL OF APPLIED PHYSICS 








ah, = or 


nm 





They produce a substantial reduction in the 
friction, they can cause smooth sliding, and they 
protect the underlying metal surfaces. In addi- 
tion, metallic films are worn off the surface by 
successive sliding over the same track in a 
manner similar to that of a lubricant film, except 
that the metallic films are worn away at a 
greater rate than the hydrocarbon films. A 
further point of striking similarity is the effect of 
melting. The transition from smooth sliding to 
stick-slips when the metallic films are melted, is 
closely analogous to the change observed when 
solid hydrocarbon films are heated through their 
melting point.® 

There are, however, several marked differences 
between metallic films and lubricant films. The 
earlier experiments showed that even on rough 
surfaces a lubricant film need only be one 
molecule thick to be effective as a boundary 
lubricant (see also Langmuir’). A metallic film 
must be appreciably thicker, of the order of 10~® 
cm if it is to be effective. A further striking and 
fundamental difference is that lubricant films 
obey Amontons’ law and metallic films do not. 
With metallic films the coefficient of friction 
decreases as the load is increased, and at high 
loads the coefficient of friction may be extremely 
low. With the indium films described in this 
paper the value of uw under heavy load (u=0.04) 
is considerably less than that observed with even 
the best boundary lubricants. This value is 
similar to that obtained on ice surfaces (see, for 
example, Bowden and Hughes’). 
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Fic. 9. Comparison of frictional behavior of: I, unlu- 
bricated steel; II, steel lubricated with mineral oil; III, 
steel lubricated with indium film of thickness 4-10~4 cm. 


7 Langmuir, J. Frank. Inst. 218, 143 (1934). 
8’ Bowden and Hughes, Proc. Roy. Soc. A172, 280 (1939). 
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Fic. 10. Microstructure of non-dendritic copper-lead alloy 


(X250). 


THE ACTION OF BEARING ALLOYS 


A series of experiments was carried out to 
investigate the part which these metallic films 
may play in the action of certain bearing alloys. 
A study was made of the friction of various com- 
binations of steel, copper, lead, silver, and 
indium and the measurements were made (a) 
with the pure metals, (6) with the metals alloyed 
with one another, and (c) with one metal de- 
posited as a thin film on the surface of the other. 
Several alloys were studied but we may here 
confine our attention mainly to copper and lead. 
Two copper-lead alloys were used. The first (see 
Fig. 10), which contained 27 percent lead was 
non-dendritic; the copper formed the con- 
tinuous phase and small droplets of the lead were 
distributed through it. The second, (see Fig. 11) 
containing 20 percent lead was dendritic; the 
lead formed the continuous phase and was dis- 
tributed between the copper dendrites. 


Frictional Behavior of Steel, Copper, and Lead 
(a) Pure Metals 


The friction of steel on copper or lead depends 
upon ‘the experimental conditions but for a 
curved steel slider on copper or on lead it was of 
the order of »~ 0.9. An examination of the steel 
surface after sliding showed the characteristic 
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Fic. 11. Microstructure of dendritic copper-lead alloy 
(X 250). 


local adhesion and welding of the softer metal on 
to the hard one. The friction of pure metals on 
one another was somewhat higher and varied 
from up~1~2, and the surfaces showed that 
some seizure and tearing had occurred. 


(6) Thin Films of Lead on Copper 


The relation between the track width and the 
friction of a curved steel contact sliding on a 
copper surface which has been coated with lead 
films of varying thickness is shown in Fig. 12. 

A film of lead 10-* cm thick (Curve II) causes 
very little reduction in the friction. As the 
thickness of the film is increased, the friction for 
any given track width decreases, and reaches 
a minimum when the film is 10-* cm thick 
(Curve V). Further increase in film thickness 
produces no further change. The frictional force 
is now governed by the shear strength of the lead 
and the width of the track and is not influenced 
by the substrate (except insofar as this may 
affect the track width). When similar measure- 
ments are made on solid lead, the points lie on 
Curve V (Fig. 12). It is clear that after the film 
has reached a thickness of about 10-* cm the 
friction is due to the interaction between the 
steel and lead. This minimum film thickness of 
10-* cm for lead films on copper is greater than 
that observed when the films are deposited on 
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steel. It is also greater than that observed with 
indium films on steel. As pointed out earlier, the 
minimum film thickness for complete ‘“‘lubrica- 
tion”’ is influenced by the hardness of the metal 


‘ substrate. With softer metal substrates, such as 


copper, which are more readily deformed, thicker 
films are necessary. It is also influenced by the 
fineness of the original surface finish and by the 
state of cleanliness of the substrate since this 
affects the strength of adhesion of the deposited 
layer. 


(c) Copper Lead Alloys 


The friction behavior of both alloys at room 
temperature was similar. Steel on the alloys gave 
smooth sliding and a coefficient of friction of 
about »=0.18. The tracks were smooth grooves 
showing little tearing. Some signs of smearing of 
extruded lead could be detected. The alloys on 
steel gave stick-slip motion, the maximum value 
of the friction being about 4.=0.3. The tracks 
showed that some metal from the upper contact 
had been welded on to the steel surface but the 
extent of this was considerably less than the 
welding observed with pure copper or lead. 

It is at once apparent that the friction of these 
alloys does not lie between the ‘values of their 
constituents. The friction of steel on copper, for 
example, is about 4.=0.9, and for steel on lead 
u=0.9. With steel on the copper-lead alloys, 
however, the friction is about ».=0.18. The alloys 
and the pure copper have approximately the 
same flow pressures at room temperature, so that 
the area of contact A is nearly the same for a 
given load. Nevertheless, the friction of the alloys 
is less than one-quarter that of their pure con- 
stituents. 

The relation between the friction F and the 
track width d for the alloy is plotted on the 
broken curve (Curve VI) in Fig. 12. It will be 
seen that these results lie very close to those 
obtained when a lead film 10~ cm thick is arti- 
ficially deposited on to copper. This at once 
suggests that the lead in the alloy is extruded 
during sliding and forms a thin “‘lubricating”’ 
film of effective thickness between 10~* cm and 
10-* cm on the copper. This is borne out by a 
microscopic examination of the track which 
showed traces of lead smeared over the surface. 
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The Effect of Wear on the Friction 


With pure metals the friction is not appre- 
ciably affected by wear. The coefficient of friction 
after the slider has passed over the same track 
20 times is similar to that observed on the first 
or second run. This is to be expected since the 
area of contact is essentially the same and the 
metallic junctions formed in each successive run 
are still of the same nature as in the first run. 

With the thin metallic films, as we have seen, 
the films are gradually worn away. As metallic 
contact occurs with the underlying substrate 
which has a much higher shear strength than the 
film materials, the friction rises. The track also 
shows the breakdown of the film and the tearing 
of the underlying metal. 

The wear results with the copper lead alloys 
are similar to those observed with lead films 
artificially deposited on to copper, and again bear 
out the view that the hard copper matrix is 
lubricated by a thin lead film. At first the ex- 
truded lead offers a reasonably effective lubricat- 
ing film. The next few runs over the same track 
may produce a more uniform smearing of the 
lead, and the friction may even decrease. The 
value of u at this stage is about 0.17. With con- 
tinued traversals over the same track, however, 
the lead film is worn away, the alloy is gradually 
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Fic. 12. Frictional behavior of lead films deposited on 
copper substrates. I, clean copper; II, lead film 10~* cm 
thick; III, lead film 10-® cm thick; IV; lead film 107-4 cm 
thick; V, lead film 10-3 cm thick and more, and solid lead; 
VI, copper-lead alloy. 
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Fic. 13. Influence of temperature on the friction of the 
non-dendritic copper-lead alloy. 


exhausted of most of the lead in the vicinity of 
the track, and the friction rises to a value of 
about 4»=0.28. This is borne out by the friction 
traces and by an examination of the track. An 
interesting confirmation of this view is to be 
found in the values obtained for the friction of 
a curved contact of bearing metal on steel. In 
this case, all the wear is concentrated on a very 
small area of the alloy which may quickly be 
depleted of lead and the friction rises rapidly to 
a value of 4.=0.28. This is almost identical with 
the value observed with steel on the bearing 
metal, after continuous running over the same 
track; i.e., it corresponds to the friction between 
steel and bearing metal surfaces when most of 
the lead between the surfaces has been exhausted. 


Effect of Temperature 


The effect of temperature on the friction of 
bulk metals and of thin metallic films has been 
discussed above. With bulk metals there is 
practically no change in friction. With thin 
metallic films the friction decreases linearly with 
temperature as the film softens. For example, 
with lead films on various metallic substrates the 
temperature coefficient of friction between 20° 
and 300°C is about y = —0.0015. 

The influence of temperature on the friction 
of the non-dendritic alloy is shown in the full 
curve (Fig. 13). The upper curved slider was of 
steel. The friction decreases in an approximately 
linear manner as the temperature is raised to 
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250°C, and the temperature coefficient y ‘is 
approximately —0.0015. The area of contact 
between the metals as determined by the track 
width is shown in the dotted curve. This shows 
little change until a temperature of 250°C is 
reached. After this the alloy softens and both the 
area of contact and the friction rise rapidly. 

Below this point the behavior resembles very 
closely that of a copper surface on which a thin 
film of lead has been deposited and it will be 
noted that the temperature coefficient y is 
approximately the same. 

With the dendritic alloy, in which the lead 
forms the continuous phase, we should expect 
the softening to occur at a lower temperature and 
to be more marked. Experiments show that this 
is the case; both the friction and area of contact 
rise rapidly above 120°C. 


The Role of Thin Films in the Action of 
Bearing Metals 


We see, therefore, that the frictional behavior 
of a copper lead alloy resembles very closely that 
of a copper surface on which a very thin film of 
lead has been deposited. The actual value of the 
friction of the alloy is the same as that of a copper 
surface on which a lead film 10~* cm thick has 
been artificially deposited. The temperature coef- 
ficient of friction is the same in each case. The 
increase in friction with wear is also very similar 
except that, in the case of the alloy there is 
evidence that the potential supply of lead is 
greater since it may continue to be expressed 
from the alloy during the sliding process. With 
a thin film of lead deposited on copper the supply 
is, of course, limited to the amount which is 
actually present on the surface at the beginning 
of sliding. 

. These results show the important part which 
may be played by thin surface films of metal in 
reducing the friction and wear of bearing alloys 
and support the view that in certain bearing 
alloys the anti-frictional and anti-seizure proper- 
ties may be due, primarily, to the spreading of 
thin films of the soft low melting constituent, 
over the surface of the harder constituent. 

Certain solids are able to provide their own 
surface films. This is the case, for example, with 
ice.* Under the frictional heating local melting 
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of the surface layers may occur and may provide 
a lubricating film, while the bulk of the solid 
remains hard. With most pure metals it is dif- 
ficult to achieve this, because of the progressive 
softening of the metal as the temperature is 
raised which leads to a corresponding increase in 
the area of contact so that the friction is little 
affected. The problem is to keep both the area 
of contact A and the shear strength s as small as 
possible. This occurs in the case of graphite and 
as is well known the friction of this is very low. 

We see, however, that the condition of small 
s and small A may be achieved artificially by 
putting a very thin film of a soft metal on to the 
surface of a hard one. This may be done by 
deposition beforehand, or it may be achieved by 
selecting an alloy of suitable composition and 
structure so that the film is spread during the 
sliding process. 

In this connection the importance of surface 
tension is apparent. The soft low melting con- 
stituent must spread readily on the harder one 
to form a thin film. This is borne out by experi- 
ment. It was shown, for example, with lead, 
indium, and mercury films on steel, copper, and 
silver that unless the soft metal readily wetted 
the harder one, it was, under severe conditions 
of sliding, relatively ineffective as a lubricant. 

Although it is desirable that these metallic 
films should adhere to and wet the substrate, the 
converse should be true about the other moving 
surface, e.g., the journal of a bearing. Provided 
the metallic film does not break down, the smaller 
the adhesion and the less the tendency of the 
metallic film to spread on to this surface the 
better. As has been pointed out in the first paper, 
one important function of a lubricant is to reduce 
the adhesion between the two moving surfaces 
and help to reduce the effective shear strengths 
of the metallic junctions. 

Apart from these surface tension and spreading 
effects we have seen that the availability of the 
soft constituent is important. The dendritic 
copper-lead alloy, for example, may under 
certain severe conditions of sliding be more 
resistant to seizure than the non-dendritic alloy. 
In the former alloy the pockets of lead are all 
interconnected so that potentially a greater 
supply of lead is available to any particular area 
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of the surface which may be in need of it. With 
the non-dendritic alloy the pockets of lead are 
isolated so that a local exhaustion of thelead 
may occur and a seizure onto the copper surface 
may take place. This difference is shown by the 
“sweating” out of the lead. When the non- 
dendritic alloy is heated, the lead appears as a 
fine mist of droplets over the surface. When the 
dendritic alloy is heated the lead collects in one 
or two large drops, by a surface tension effect. 
The dendritic alloy thus has a three-dimensional 
supply of lead to the surface while with the non- 
dendritic alloy it is restricted to the surface layers. 

In this paper the experiments and discussion 
have deliberately been confined to a considera- 
tion of the frictional and seizure properties of 
the metals. It is clear that the performance of a 
bearing alloy must also be dependent on its bulk 
properties, such as mechanical strength, thermal 


conductivity, fatigue resistance, etc., and unless 
these are satisfactory the bearing will not func- 
tion. The influence of the mechanical and bulk 
properties of bearing metals has been studied by 
numerous authors. A recent discussion of some 
of them has been given by Gillet, Russell, and 
Dayton.?® 

In this paper the discussion has been confined 
to unlubricated metals. Even in the presence of 
a lubricant, metallic contact between the moving 
surfaces may readily occur so that the basic fric- 
tional properties of the metals themselves are of 
primary importance. Nevertheless, the whole 
behavior is profoundly modified by the presence 
of lubricants and by the extent to which the 
surface layers of the various metals are capable 
of adsorbing the lubricant films. 





® Gillet, Russell, and Dayton, Metals and Alloys 12, 
274 ff. (1940). 





New Formulas for Calculating Short-Circuit Stresses in Bus Supports for 
Rectangular Tubular Conductors 


THOMAS JAMES HIGGINS 
Illinois Institute of Technology, Chicago, Illinois 


(Received December 7, 1942) 


New and simple formulas are derived for calculating short-circuit stress in bus supports 
for rectangular tubular conductors. A numerical example illustrates their use and advantages. 


HE mutual electromagnetic force exerted 

between two long isolated nonpermeable 
conductors, spaced D inches between their 
parallel axes and carrying currents of J and I’ 
amperes, is commonly calculated from the for- 
mula 


f=(5.4 RII'D) 10-7 pound per foot of bus length. 


Herein k is a parameter, commonly termed the 
electromagnetic space factor, of magnitude deter- 
mined by the geometry and relative position of 
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the two conductor cross sections. As each of the 
several methods available for calculating the 
mechanical force acting on bus bar supports 
during short-circuit requires a priort knowledge 
of the pertinent values of k, it is evident that of 
vital interest to the designer of single phase or 
polyphase bus structures are the values of k for 
the type of bus most commonly used: that 
comprised of a group of identical conductors 
having parallel, coplanar axes and parallel-sided, 
full or hollow rectangular cross sections. 
Formulas and curves for conductors of full 
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rectangular cross section were presented by H. B. 
Dwight! a quarter of a century ago. Only re- 
cently, however, have formulas been presented 
for conductors of hollow rectangular cross section 

a type being used in increasing amount to 
carry the very heavy current loads encountered 
in new defense plants that house large units for 
welding, induction heating, electric smelting, 
electrolytic processing, or kindred operations. 

Elsewhere the has shown that for 
identical conductor cross sections located as in 
Fig. 1 the electromagnetic space factor k is given 
by 


writer? 


k=[D?/(rs—RS)?][M(r, s) +M(R, S) 
+2M(R+4t, t)+2M(t, S+2) 
—2M(R+t, S+t)—2M(t, t)], (1) 


wherein, typically, M(r, s) is defined by 


M(r, s)/r?=2m tan—! m—log (1+?) 


9 9 








—143m? 
+, (2) 


n*m? n' 
$4 
6(1+m?*) 90 (1+-m?)8 
and m=s/D,n=r/D. 
A form of M(r,s) especially convenient to 
rapid numerical computation is obtained by 
expanding the terms of the right-hand member 


1H. B. Dwight, Elec. World 70, 522 (1917). 
2 T. J. Higgins, Trans. A.I.E.E. 61, 578 (1942). 
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of (2) in powers of m: 


u 


n> n‘ 
M(r, s)/ =m (1 oie -) 
6 15 


1 nm 14n* 
+mt(—+"4— Jae. (3) 


Commonly, the geometry of a given problem is 
such that in (2), and hence in (3), the term in n* 
is negligible. Often, the terms in both n? and n‘ 
are negligible. 


Example No. 1 


Illustrative of foregoing remarks, if k is sought 
for two 2.5 by 2.5 inches square tubular con- 
ductors, 0.5 inches thick, spaced 10 inches 
between parallel coplanar axes, (1) yields 


k=(100/16)[M(2.5, 2.5) +M(1.5, 1.5) 
+2M(2, 0.5)+2M(0.5, 2) 
—2M(2, 2)—2M(0.5,0.5)] (4) 
= (100/16) [0.39049+0.05062+0.02012 
+ 0.01988 — 0.31992 — 0.00126] 
= 7.999 + 


wherein, typically, we obtain from (3)—the 
terms in m‘ are neglected, for in (2) the term in 
n' is negligible— 


M(2.5, 2.5) =6.25(0.063151 — 0.000692 
+0.000019) =0.39049. (5) 


This example illustrates both the advantages 
and disadvantages of using (1), (2), and (3) in 
the form given. As evidenced by (5), the right- 
hand member of (3) converges rapidly, whence 
only a few terms need be calculated. On the 
other hand, these terms must be calculated 
correct to many decimal places; for, as evidenced 
by (4), the value of the right-hand member of 
(1) is determined primarily as the sum of terms 
of small magnitude and of a difference of the 
same order of magnitude. Accordingly, as the 
number and type of computations required by 
(1) are most conducive to numerical error, and 
as k must be determined with accuracy if it is to 
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possess significance whatsoever, much simpler 
formulas that yield k with but little computation 
are obviously desirable. Such formulas are pre- 
sented in this paper. 


1. DERIVATION OF THE NEW FORMULAS 
We obtain from (3) 


rs? s2 
M(r, s)= ( al —) 
D>? 6D* 15D8 


rs? 


righ 


r2s4 posh 
-(—+—+—) 
6D* 6D® 6D* 


r2s§ — rts6 = 1 456 
+ (—- + ;) . (6) 
AI5SD® 6D® 45p1 





Substituting appropriate forms of (6) in (1), 
discarding the term containing D® as negligible, 
and arranging the remaining terms as a series in 
powers of D yields 


1 
k=(rs—RS)*} > r’s?+-—[D (rs? —r’s*) ] 
6D? 


1 
+—[) (2r's? — Sr's*+-r°s*) ] 
30D' 


1 
+—[> (rst+r's*)]}, (7) 
6D* 


wherein the summation signs indicate manipu- 
lations typified by 


7s? =(7sk+ RS +2(R+1)°24+2f(S+t)? 
—2(R+4#)?(S+i)?—2f] (8) 
=42(N+M)?, 


M and N defined by M=S+t and N=R+41. 
Similarly, 


Dr's? =42[N2(N2+4NM+43M2) 
+£(M?2+4MN+3N?)], (9) 


Dr's? = 2PLN4(2N2?+12NM+15M?) 
+5?N?(3N?+8NM+15M?) 
+#(2M?+12MN+15N?) ], (10) 
Drs! =4P[M2(M?+4MN+3N?) 
+?(N?+4NM+3M?) |, 
Dr's! = 4P(N?M?(3N?+8NM+3M2) 
+?(N'+8N?M+18N?M?+8NM*+M‘) 
+#(3M?+8MN+3N?")], (12) 


(11) 
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Dr's! = 22[3N'M2(2N2+8NM+5M?) 
+2N2(2N'+24N°M+90N?2M? 
+80NM?+15M*)+1(2M! 
+24M*N+90M?2N?+80M N?+15N') 

+31°(2M2+8MN-+5N?), (13) 
Dr's* = 2e[ M4(2M?2+12MN-+15N?) 
+52M*?(3M?+8MN-+3N?) 
+1(2N2+12MN+15M2)], (14) 

Dr's = 2°[3M*N2(2M?2+8MN+5N?2) 
+PM2(2M!+24M*N+90M2N? 
+80MN?+15N*)+14(2N' 
+24N*M+90N?M?2+80NM?+15M!) 

+31°(2N?+8NM+5M2)]. (15) 


Substituting (8) to (15) inclusive in (7) and col- 
lecting terms yields 


k=(1/(N+M)*]{(N+M)? 
+ (N?— M?)(N?+ M?+4NM-+2?) /6D? 
+[(N®+6N'M —20N°M? 
+6NM*®+ M*)+50(N*—6N?M?+ M‘*) 
+t1(N?—8NM+ M?) |/15D* 
— (N?— M*)[6(N?+4NM-+ M?) M?N? 
+2#(N‘+24N°M+76N?M 
+24N M*+ M*)+#(13N?+56NM 
+13M*)+7t®]/12D*}. (16) 
The hollow rectangular conductors currently 
used are, for the most part, square. For this prac- 
tical important case M= JN, and (16) reduces to 
the absurdly simple formula 


k=1—(3M'+10M*? +3) /30D* 
=1—(3M?+?)(M?+3)/30D‘. (17) 


2. AN ILLUSTRATIVE PROBLEM 
Example No. 2 


Utilizing the same data as in example 1 
whence M=2 inches, t=0.5 inches, and D=10 
inches—we obtain from (17) 4 


k=1—(48+10+0.1875) /300,000 
=1—0,000019=0.9999+. (18) 
Contrast of the computation of example 2 with 
that of example 1 emphasizes, more strongly 
than can words, the advantages stemming from 
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the use of (16) or (17) rather than (1), (2), and 
(3). 

Corroborative of the foregoing analysis: if 
R=S=0 and r=s, we have conductors of full 
square cross section and (1) reduces to 


k=(D*/r')M(r, 1); 


whence substitution of M(r,r) from (3) 


after 


discarding the term in r'°—yields 
k=(1—(M*/30) ]. 


But this is precisely what we obtain from (17) 
on substitution of M=t=r/2. 

Finally, we note that all expressions for k are 
dimensionless; whence the linear dimensions can 
be expressed in units chosen at will. 
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